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This review examines intermittent target search strategies, which combine phases of slow motion, 
allowing the searcher to detect the target, and phases of fast motion during which targets cannot 
be detected. We first show that intermittent search strategies are actually widely observed at 
various scales. At the macroscopic scale, this is for example the case of animals looking for food ; 
at the microscopic scale, intermittent transport patterns are involved in reaction pathway of DNA 
binding proteins as well as in intracellular transport. Second, we introduce generic stochastic 
models, which show that intermittent strategies are efficient strategies, which enable to minimize 
the search time. This suggests that the intrinsic efficiency of intermittent search strategies could 
justify their frequent observation in nature. Last, beyond these modeling aspects, we propose 
that intermittent strategies could be used also in a broader context to design and accelerate 
search processes. 
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I. INTRODUCTION 

A. General scope and outline 

What is the best strategy to find a missing object ? Anyone who has ever lost his keys has already been facing 
this problem. This every day life situation is a prototypical example of search problem, which under its simplest form 
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involves a searcher - either a person, an animal, or any kind of organism or particle - in general able to move across 
the search domain, and one or several targets. Even if very schematic, the search problem as stated turns out to be 
a quite universal question, which pops up at different scales and in various fields, and has generated an increasing 
number of works in recent years, notably in the physics community. 

Theoretical studies of search strategies can be traced back to World War II, during which the US navy tried to 
most efficiently hunt for submarines and developed rationalized search procedures ( Champagne et a"ll . l2003tlShlesingei1 . 
I2009D. Similar search a lgorithms have then been developed and utilized in the context of castaway rescue operations 
( Frost and Stonel . l200lh . or even for the recovery of an atomic bomb lost in the Mediterranean sea near Palomares 
i n 1966. One can cite for in stance the rescue of the Scorpion, a nuclear submarine lost near the Azores in 1968 
( Richardson and Stond . Il97l[ ). At the macroscopic scale, other imp ortant and wi dely studied ex a mples of search 
processes conc e rn an i mals searching for mate, food o r shelter (iBelj [19 91; Bcnic hou et al 1 l2006bl : ICharnovl . Il976l : 
lEdwards etldl l2007l : lO'Brien et ail . Il990t IShlesingeij . 120061 : IViswanathan et al.i . |1999[ ) which will be discussed in 
more details in this review. One can even mention pre-historic migrations, which, apart from classical archaeological 
literature, ha ve also been studied as a search problem, in which human groups search for new profitable territories 
(|Floresl . [2007h . At the microscopic scale, search processes naturally pop up in the context of chemical reactions, for 
which the encounter of reactive molecules - or in other words the fact that one searcher molecule finds a reactive 
target site - is a required first step. One sho uld obviously menti o n the theory of diffusion-controlled reactions, initiated 
many years ago by the cele brated theory of 'von Smoluchowskil ( 1917t ) and developed by innumerable researchers (see 
for instance the review bv iHanggi et al". (19.90.)). More recently, this field has regained interest in the context of 
biochemical reactions in cells, where the sometimes very small number of reactive molecules makes this first step of 
search for a reaction partner crucial for the kinetics. One can think for instance of reactions involved in genomic 
transcription, a representative example of which being the search for spec i fic DNA sequences by tr anscription factors 
(|Berg et al.l , ll98ll:lBonnet et all l2008l : iGorman and Greend . 120081: iMirnvl . l2008l : l^nHiiiel, [2003) . 

In all these examples, the time needed to discover a target is a limiting quantity, and consequently minimizing this 
search time often appears as essential. In order to get an intuition of what could be an efficient search strategy on 
general grounds, let us go back to the everyday- life example mentioned above (see figure [Ij. We consider a searcher 
who lost a tiny object - let us say a key - in a large sandy beach, where the key is so small that it cannot be detected 
if the searcher passes by too fast. In addition, we assume that the searcher has no prior information on the position 
of the key, except that the key is in a bounded domain (the beach) . What is then the best strategy for the searcher to 
find the key as fast as possible? A first strategy consists in a slow and careful exploration (to make sure that the key 
will be detected upon encounter) of the sand all along the beach. In the case of a very large beach, the search time can 
then be very long. An alternative strategy one can think of consists in interrupting the slow and careful exploration 
of the sand by mere displacement phases, during which the searcher relocates on the beach very fast, but without 
even trying to detect the key (typically the searcher "runs"). We will call hereafter intermittent search strategies 
such processes that combine two distinct phases : a phase of slow displacement which enables target detection, and 
a phase of faster moti on during which th e target cannot be detected (note that the word intermittent has also been 
used more recently by iBartumeu^ ( 20091 ) with another definition) . 

The efficiency of such intermittent strategies results from a trade-off between speed and detection and can be 
qualitatively discussed. Intuitively, the advantage of the fast relocation phases for the searcher is to reach unvisited 
regions. The drawback is however that during these phases time is consumed without any chance of detecting the 
target. The net efficiency of this strategy is therefore not trivial, and these last years, many works have focused on the 
following questions: (i) Can phases of fast motion which disable detection make the global search more efficient? (ii) 
If so, is there an optimal way for the searcher to share the time between the two phases ? (iii) Are these intermittent 
search patterns relevant to the description of real situations ? 

The goal of this article is to review these works while bringing explicit answers to these questions. More pre- 
cisely, we will first show that intermittent transport patterns are actually widely observed at various scales. At the 
macroscopic scale, this is for example the case of foraging animals (see section |lT| ; at the microscopic scale, we will 
show that intermittent transport patterns are involved in reaction pathway of DNA binding proteins as well as in 
intracellular transport (see section llll| . Second, we will show on generic stochastic models that intermittent strategies 
are efficient strategies which enable to minimize the search time (see sections IIIIIIIIIV]) , and therefore suggest that 
this efficiency could justify their frequent observation in nature. Last, beyond these modeling aspects, we will propose 
that intermittent strategies could be used also in a broader context to design and accelerate search processes. 



B. General framework and first definitions 



The search problem can take multiple forms ( da Luz et al.l l2009l ): in this section we define more precisely the 
framework of this review - namely the random intermittent search strategies - and introduce the main hypothesis 



Figure 1 Illustration of intermittent reaction paths by an every-day life example of search problem. The searcher looks for a 
target. The searcher alternates fast relocation phases, which are not reactive as they do not allow for target detection, and 
slow reactive phases which permit target detection. 



which will be made. 



1. Searching with or without cues 

Although in essence in a search problem the target location is unknown and cannot be found from a rapid inspection 
of the search domain, in practical cases there are often cues which restrict th e terr i tory t o explore, or give indications 
on how to explore it. We can quote the very classical exampl e of chemotaxis (lBere , l2003) . which keeps arising i nteres t 



in the biological a nd ph ysical communities (see for exa mple lKafri and Da Silveiral l|2008D ; iTailleur and Gated (|2008l ): 
lYuzbasvan et al.l (j2003l )). Bacteria like E.coli swim with a succession of "runs" (approximately straight moves) and 
"tumbles" (random changes of direction). When they sense a gradient of chemical concentration, they swim up or 
down the gradient by adjusting their tumbling rate : when the environment is becoming more favorable, they tumble 
less, whereas they tumble more when their environment is degrading. This behavior results in a bias towards the most 
favorable locations of high concentration of chemoattractant which can be as varied as salts, glucose, amino-acids, 
oxygen, etc... More recently it has bee n shown that a similar behavior can also be t riggered by other kinds of external 
signals su ch as temperature gra dients ( Maeda et al.l , Il976l : ISalman and Libchabeil l2007l : I Salman et al 1 120061 ) or light 



intensity ( Sprenger et al.l , [l993f l. 



Chemotactic search requires a well defined gradient of chemoattractant, and is therefore applicable only when the 
concentration of cues is sufficient. On the contrary, at low concentrations cues can be sparse, or even discrete signals 
which do not allow for a gradient based strategy. It is for example the case of animals sensing odors in air or water 
where the mixing in the potentially turbulent flo w breaks up the chemical signal into random and disconnected patches 
of high concentration. IVergassola et al.l ( 2007f) proposed a search algorithm, which they called 'infotaxis', designed 



to work in this case of sparse and fluctuating cues. This algorithm, based on a maximization of the expected rate of 
information gain pr oduces trajectories such a s 'zigz agging' and 'casting' paths which are similar to those observed in 
the flight of moths ( Balkovskv and Shraimanl . |2002| ) . 



In this review we focus on the extreme case where no cue is present that could lead the searcher to the target. This 
assumption applies to targets which can be detected only if the searcher is within a given detection radius a which 
is much smaller than the typical extension of the search domain. In particular this assumption clearly covers the 
case of search problems at the scale of chemical reactions, and more generally the case of searchers whose motion is 
independent of any exterior cue that could be emitted by the target. 



2. Systematic vs random strategies 

Whatever the scale, the behavior of a searcher relies strongly on his ability, or incapability, to keep memories of his 
past explorations. Dependin g on the sea rcher and on the space to explore, such kind of spatial memory can play a 
more or less important role ( Moreau et al. , 2009ah . In an extreme case the searcher, for instance human or animal. 



can have a mental map of the exploration space and can thus perform a systematic search. Figure 2 presents several 
system atic patterns : lawn-mower, expanding square, spiral (for more patterns, see for example IChampagne et ahl 
( 20031 )). These type of search have been extensively studied, in particular for designing efficient search operated by 



humans (Dobbie, 1968; Stone, 1989). 

In the opposite case where the searcher has low - or no - spatial memory abilities the search trajectories can be 
qualified as random, and the theory of stochastic processes provides powerful tools for their quantitative analysis 
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Spiral Expanding square Lawn-mower 



Figure 2 Examples of patterns for systematic exploration of space 



(for a short reminder on random walks, see Appendix |A|. This is obviously the case for "molecular" searchers at 
the microscopic scale which are subject to thermal Brownian motion, but also at larger scales of animals with low 
cognition skills. This review is mainly focused on random search problems, and effects of spatial memory will be 
briefly discussed in the last section. 

Note that we will keep using the word "strategy" for animals with low cognitive abilities and even for molecules, 
although such searchers are not able to design strategies themselves since of course their dynamics are simply governed 
by the laws of physics. In the context of proteins searching for targets on DNA, we mean that the search time depends 
on parameters such as the ionic strength or the protein/DNA affinity which, if varied, can lead to a minimization of 
the search time. In the case where the search kinetics is a limiting constraint, such good or even optimal values of 
these parameters might have been selected in the course of evolution. We call strategy this very fact that physical 
parameters can be tuned (implicitly by evolution) to optimize a biological function. Note however that the real 
optimization problem depends on many parameters and constraints. The models studied in this review are restricted 
to kinetic constraints, which can be dominant both at the microscopic and macroscopic scale, as discussed in sections 
II and III. This key assumption will be used throughout this review. 



3. Framework of this review 

To summarize, in this review we shall focus on intermittent search strategies for targets which emit no cue. The 
searchers will be assumed to have no (or low) memory skills, resulting in their trajectories being intermittent random 
walks. Depending on the example to be treated, different quantities can be used to assess the efficiency of search 
strategies, such as the energy necessary for reaching the first prey, the number of preys collected in a given time, or 
the time for encountering the first prey. In this review we will discuss the efficiency of search strategies uniquely from 
a kinetic point of view. We shall mainly consider the mean first passage time to a target as a quantitative measure of 
the search efficiency, and study the minimization of this quantity. Note that the full distribution of the first-passage 
time is a priori needed to quantify the search kinetics on all time scales. However, in most of the situations that 
we will consider in this review, in can be checked numerically that the distributions of the search time can be well 
approximated by an exponential, which means that the kinetics is fully characterized by the mean first-passage time. 



II. INTERMITTENT SEARCH STRATEGIES AT THE MACROSCOPIC SCALE 



Searching for a randomly l ocated obje ct is one of the most frequent tasks of living organisms, be it for obtaining food, 
a sexual partner or a shelter ((Bdl .[ 199ll) . In these examples, the search time is generally a limiting factor which has to 
be optimized for the survival of the species. The question of determining the ef ficiency of a search behavior is thus a 
crucial problem of behavioral ecol ogy, which has ins pired numerous e xperimental (iBell ^ 1991 [ Kramer and McLaughlin , 



2001 



2008 



lO'Brien et al.L Il990l. Il989l) and theoretical iBenichou erall l2005bl l2006bt iBover et al.l l2006l : iLomholt et al. 
IViswanathan et al.Ul996l . Il999l ) works . In this context. Levy walk strategies have been proved to play a crucial 



role in such optimization problems. In this section, we first remind why these Levy walks are advantageous with 
respect to simple random walks when searching randomly, as first put forward by IShlesinger and Klafted (I1986D . We 
also remind the pioneering model of Viswanathan et al., which has played a major role in the development of ideas on 
random search strategies. On the other hand, we show how intermittent strategies are naturally involved as soon as 
hidden targets are considered and define a basic model relying on intermittent strategies, introduced to account for 
the search behavior of "saltatory" animals. This one-dimensional model is then extended to a bi-dimensional model. 
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which is shown to be a minimal model optimizing the search time. Last, we discuss the relationships between these 
two main classes of search strategies - Levy and intermittent - and come back on the famous "albatross story". 



A. The Levy strategies 

1. The advantage of Levy walks with respect to simple random walks 

The ballistic pha s es intersper sed with t urns of animals traject o ries h ave often been interpreted as Levy walks 
(jViswanathan et all Il999l 120081) . Actua lly, [Shlesinger and Klafted (|1986[ ) were the first to report that, due to their 



weak oversampling properties (see figure lILA.ip , Levy walks could be an efficient way to explore space and could be 
used to model in particular trajectories of foraging animals. As a matter of fact, the mean number of distinct sites 
visited in n steps - which is a measure of the territory explored - is known to behave for a standard random walk 
like n in dimension d> 2, and like n'^/^ if d < 2. This is less efficient in low dimension than for a Levy walk, of jump 
probability in dimension d of the form p{r) oc r~^~^ (where (3 is the index of the walk), for which the mean number 
of distinct sites visited in n steps behaves like n, as soon as /3 < d. From the point of view of the extension of the 
territory explored after a given number of steps, the advantage of Levy walk patterns over standard random walks is 
thus clear, and this effect turns out to be as much marked as the number of searchers involved in the process is high 
( Viswanathan et ari , ll996[ ). 



Figure 3 Example of Levy walks, with /i — 1.5 (green, not present on the zoom), /i = 2 (red), jj. — 3 (blue). The total path 
length is the same for the 3 examples. 



2. Optimizing the encounter rate with Levy walks : how and when? 

These observations led IViswanathan et aL 1 (ll999^ to propose the following Levy search model, in the presence of 
fixed targets randomly and sparsely distributed : they consider a searcher performing ballistic step at constant speed, 
and detecting targets closer than r„. A target is found when the searcher encounters it for the first time. The step 
lengths are drawn from a Levy distribution p(l) cx , with 1 < ^ < 3. For /i < 1, the probability distribution is not 
defined. For 1 < /i < 2, the distribution has no mean and no variance. For 2 < < 3, the distribution has a mean 
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but no variance. For /i > 3, the distribution has both a mean and a variance, thus it obeys the central hmit theorem : 
after enough steps, the probability distribution of the difference between the starting point and the last position is a 
G aussian, as if it were diffus ion, with the mean square distance scaling linearly with time. 
IViswanathan et al. ] (Il999f ) are interested in the mean number of targets detected after a large observation time t. 



More precisely, they ask the following question : is it possible to optimize this number with respect to the exponent 
fi characterizing the motion of the searcher ? To answer, they actually consider two different types of targets, which 
lead to two very different optimal strategies. 

• In the first case of what they call "revisitable targets" - meaning that, as soon as detected, a target reappears 
at the same location - they rely on a mean-field approximation of the problem and find that the encounter rate 
is optimized for a Levy exponent ^ ~ 2. 

• In the second case of "non revisitable targets" (or destructive search) where e ach target can be found onl y once, 
or in the case of a single available target, the optimal strategy proposed in IViswanathan et al.l ( 1999^ is not 
anymore of Levy type, but reduces to a simple linear ballistic motion. 

Several extensions of th is pioneering model have been proposed. 

iBartumeus et al.l ( 20021 ) studied the case of non-revisitable mov ing targets. They s howed that a Levy strategy with 



/i = 2 is often better than a "Brownian" one > 3). However. iJames et al.l ( 20081 ) extended the study to ballistic 
motion, that outperformed these Levy strategies. 



An interm e diate situation between revisitable and non revisitable targets has been studied bv lRaposo et al.l ( 2003f l: 



ISantos et al.l ( 2004 l. In these works, the immobile target is destroyed upon encounter, but regenerates after a time 



T at the same place (for example a plant bearing new fruits after previou s fruit s have been eaten). Two regimes 
are found. When r is large (> Tc, a critical time evaluated in lRaposo et al.l ( 20031 ): ISantos et al.l ( 2004 )). the simple 



ballistic motion remains the best strategy. When r < Tc, the best fi is between 1 and 2. However, one could argue 
that in this regime the simple strategy where the searcher does not move and waits for the renewal of the target 
outp erforms searching for an hyp othetical other target. 



In IBartumeus and LevinI ( 20081 ). the targets are in patches (such as fish schools) or Levy distributed. Even if the 



targets are destroyed upon encounter, finding a target means that the presence of other targets in the vicinity is likely, 
which is close to the case of revisitable targets. Hence, like for revisitable targets, the optimum is achieved for a Levy 

distr ibution, with ii 2. 

In iRevnolds and Bartumeuj ( 2009t ). the optimum for destructive targets is /i — > I except in two cases (where 



1 < < 2). On the one hand the optimum is not ballistic when the searcher can fail in capturing a detected 
target. On the other hand, for targets destroyed upon encounter, and for the very specific one-dimensional case, as 
the measure of efficiency is the number of targets captured during a long time, the searcher is after some time in a 
situation with a target close on one side, but the next target on the other side very far away : a pure ballistic motion 
is not favored because it can take the wrong direction. 

Finally, in the case of revisitable targets and the related cases (regenerating targets, patches, failed capture), the 
Levy strategy fi = 2 emerges as a compromise between trajectories returning to one and the same ever target zone, 
and straight ballistic motion which is indeed the best way to explore space. Note however that, as stated above, in 
this case the strategy which consists simply in waiting for target renewal performs even better. In the case of non 
revisitable targets - the generic situation considered hereafter - the best strategy for the searcher is a mere ballistic 
motion without reorientations. 

In all these Levy walks models, the searcher is assumed to be able to detect targets all along its trajectory. 
Qualitatively, it corresponds to the case of targets "not too difficult" to detect. However, as it was the case with the 
example of the lost small key given in introduction, it is evident that in some situations the velocity degrades the 
perception. What happens if the targets are really "hidden", that is to say more precisely, if searching and moving 
are incompatible ? In recent years, many works have been devoted to answering this question. Most of them rely 
on the following simple two-state model, historically introduced to account for the search behavior of the "saltatory 
animals". 



B. A basic model of intermittence 

1. Observations: tfie case of "saltatory animals" 

Anyone who has ever lost his keys knows that instinctively we often adopt an intermittent behavior combining 
local scanning phases and relocating phases. Indeed, numerous studies of foraging behavior of a broad range of 
animal species show that such an intermittent behavior is commonly observed and that the durati ons of search and 
displacement phases vary widely ([Belt Il99lt iKramer and McLaughiinl . l200ll: lO' Brien et al. I Il990f) . The spectrum, 
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which goes from cruise strategy (for large fishes that swim continuously such as tuna), to ambush or sit-and-wait 
search, where the forager remains stationary for long periods (such as rattlesnake), has remained uninterpreted for a 
long time. As explained in the introduction, the interest of thi s type of intermittent strategies, often referred to as 
"saltatory" (jKramer and McLaughlii] . l200ll : IO'Brien et al.l . [l990l) in the context of foraging animals, can be understood 
intuitively when the targets are "difficult" to detect and sparsely distributed, as it is the case for many foragers (such 
as ground foraging birds, lizards, planktivorous fish...^): since a fast m ovement is known to significantly degrade 
perception abilities ( Kramer and McLaughlinl 1200X1 : lO'Brien et al. 1. 119901) . the forager must search slowly. Then, it 
has to relocate as fast as possible in order to explore a previously unscanned space, and sear ch slowly again. 

Even though numerous models based on optim ization of the net energy gain ( Anderson et al.l . Il997t 
iKnoDpien and Reddingiu"E Il985t lO'Brien et al.l . Il989t ) predict an optimal strategy for foragers, the large number 
of unknown parameters used to model the complexity of the energetic constraint, renders a quantitative comparison 
with experimental data difficult. In the model presented in this section, the search time is assumed to be the relevant 
quantity optimized by the forager in order to obtain a sufficient daily amount of food and to precede other competing 
foragers. The energy cost is treated only as an external constraint that sets the maximal speed of the animal. As 
explained in the next sections, this purely kinetic mod el of target search captures the essential features of saltatory 
search behavior observed for foragers in experiments ( Kramer and McLaughli"nl . l200l[ ). when the predator has no 
information about the prey location. 



2. The Model 



Move . 
state 2 

Search 
state 1 




Figure 4 Basic model for intermittent search 



The central point of this schematic model ( Benichou et al.l , l2005b( ) is that it relies on the explicit description of 



searching trajectories as intermittent. In the following it is assumed that the searcher displays alternatively two 
distinct attitudes (see figure 2]) : 

• a scanning phase, named phase 1, during which the sensory organs of the searcher explore its immediate vicinity. 
This phase is modeled as a "slow" diffusive movement (a continuous random walk with diffusion coefficient D). 
The target is found when this movement reaches the target location for the first time. As focusing and processing 
the information received by sensory organs require a minimum time, the scan phase cannot be too short, which 
implies a minimal mean time spent in this phase r™™. 



a motion phase, named phase 2, during which the searcher moves "fast" and is unable to detect targets. These 
relocating moves are characterized by a ballistic motion (at constant velocity V). In the case of an imals, there 
is usually correlations in the angles between two successive ballistic phases ( O'Brien et al.1 , Il990[ l. We limit 



Note that there are counter-examples such as birds of prey which can detect targets even at large velocities 
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ourselves here to the case of high correlations, which allow us to consider an effective 1-dimensional problem for 
both phases, with phase 2 always in the same direction. 

Next, it is assumed that the searcher randomly switches from phase 1 (respectively, 2) to phase 2 (respectively, 1) with 
a fixed rate per unit time Ai (respectively. A2); that is with no tempo ral memory. It leads to exponent i ally d i stributed 
phase durations, in agreement with numerous experimental studies (jFuiiwara et al. 1 120021 : [mil et al. I I2000t iLi et iZl 
I2OO8I : [Pierce-Shimonura et al.l [TqQS ) . the mean duration of phase i being = l/A^. Last, the preys are assumed to 
be immobile (see section lV.DI for a discussion of moving vs immobile targets). 



3. Equations 

We now evaluate the average time needed to find a target. The chosen geometry is a single target in a; = on 
a segment of size L with periodic boundary conditions. This geometry is equivalent to the case of regularly spaced 
targets or to the case of one target centered in a finite domain with reflective boundaries. L is thus the typical distance 
between targets, or the size of the search domain. The instantaneous state of the searcher can be described by its 
position X on the segment and by an index i which specifies its motion: 1 corresponds to the slow detection phase, and 
2 to the ballistic non-reactive phase. The survival probability pi(t, x) that, when the searcher starts at time t ~ Q from 
X and in state i, the t arget has not yet been found at time t, is known to satisfy the backward Chapman-Kolmogorov 
differential equations ( Gardineil . Il996l : iRedneij . l200lh : 



Knowing that ti{x), the mean first passage time at the target, starting from x in phase i, is given by 

U{x)^- [°°t^^%^dt^ r p,{t,x)dt, (3) 
it is easily found from Eqs.(IT]), © to satisfy 



dt JO 



D^ + -{t,{x)-h{x)) = -l (4) 
ax Ti 



-V^ + -ihix)~h{x))^-l (5) 

ax T2 

These differential equations have to be completed by boundary conditions. As we have periodic boundary conditions 
and as the target at x = can be found only in state 1, we get ti{0) = ti{L) = 0, ^2(0) = t2{L). 



4. Results 



The average search time {t) is defined as the average of ti{x) over the initial position x of the searcher, which is 
uniformly distri buted over the segment [0, L], as the searcher initially does not know the target's location. It is found 
to be given by ( Benichou et al.l , [2005b( ) : 



(t> = (T2+Tl) 



/ L (e"+'' - 1) yrrif + (1 + 2r) (e'^ - e") 1 \ 
1^2" ^/l + Ar {cP - 1) (c" - l)T2y r j ' 



with 



(6) 



(7) 



"=tr"r2V^ 



(8) 
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In the limit of L » Vt2, VDt\, Dti/{Vt2), this expression simphfies : 

2t2V ^JDTl^/IT^^TT^2^' 

Note that thanks to intermittence, {t) oc L, whereas for diffusion alone the mean detection time is tdiff = i^/(12£')). 
Intermittence is thus favorable (meaning that the gain, defined as tdiff/ (t) is grater than 1), at least for L large enough. 

Intermittence is favorable and the strategy can even be optimized. The mean search time is minimized for r^*^* = 
t{"™, and Tj^* satisfying the relation : 

^3 + 6^-84=0, (11) 

where r = D/V"^ is an extra characteristic time, depending on the searcher's characteristics. This minimum takes a 
simple form in two different regimes. 



• If Ti 3> T, the minimum of the search time is for ri = t["™, and : 

(12) 



opt _ I -Ji 'l 



In this regime, denoted by S for "searching", one has ti > T2 : the searcher spends more time scanning than 
moving . 

• If Ti <C T, the minimum of the search time is for ti = 7-™™^ and : 

opt 



(13) 

In this regime, denoted by "M" for moving, one has ti < T2, which means that the searcher spends more time 
moving than scanning. 



5. Comparison with experimental data 

The se results have been compared to experimental data from lO'Brien et all (1993) and iKramer and McLaug hlinl 
(|200lD . who provide the average duration of detection and ballistic phases, characterizing the saltatory behavior of 
18 different species, as various as planktivorous fish, ground foraging birds, or lizards. The op timal strategy obtained 
above is shown to account reasonably well for these data (see fig. [5]and lB"enichou et al. ] (|2005bll3 ldl) for further details). 

These results show that the saltatory patterns observed are a way to optimize the search, and that it is probably a 
reason why this type of patterns is observed so often, as it could have been favored by natural selection. 



C. Two-dimensional intermittent search processes: An alternative to Levy strategies 

1. Motivation 

The model of intermittent search presented previously was one-dimensional, with ballistic phases infinitely cor- 
related, in the sense that t he direction taken is always the same. Here we present a model of intermittent search 
strategies in dimension two (jBenichou et al ]. l2006bL[2007h which encompasses a much broader field of applications, in 



particular for animal or human searchers. It is shown that bidimensional intermittent search strategies do optimize 
the search time for non revisitable targets, i.e. targets that are destroyed upon discovery (see section III.A.2[) . The 
optimal way to share the time between the phases of non reactive displacement and of reactive search is explicitly 
determined. Technically, this approach relies on an approximate analytical solution based on a decoupling hypothesis, 
which proves to reproduce quantitatively numerical simulations over a wide range of parameters. 
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phase 1 




phase 1 

,--f 



'"\ phase 2 



Static mode. The slow reactive phase is static and 
detection takes place with finite rate k. 



phase 2 




Diffusive mode. The slow reactive phase is diffusive and 
detection is infinitely efficient. 



Figure 6 Two models of intermittent search: The searcher alternates slow reactive phases (regime 1) of mean duration ri, and 
fast non reactive ballistic phases (regime 2) of mean duration T2 



2. The Model 



Following the previous model, we consider a two state searcher (see figure ^ of position r that performs slow 
reactive phases (denoted 1), randomly interrupted by fast relocating ballistic flights of constant velocity V and 
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random direction (phases 2). We assume the duration of each phase i to be exponentially distributed with mean r;. 
As fast motion usually strongly degrades perception abilities, we consider again that the searcher is able to find a 
target only during reactive phases 1. The detection phase involves complex biological processes that we do not aim 
at modeling accurately here. However, we put forward here two modes of detection. The first one, referred to in the 
following as the "diffusive mode", corresponds to a diffusive modeling (w ith diffusion coefhcient D) of the search phase 
like i n the previous model, in agreement with observations for vision (iHuevl . [ 1968( ). tactile sense or olfaction 
Il99ll ) . The detection is assumed to be infinitely efficient in this mode : a target is found as soon as the searcher-target 
distance is smaller than the reaction radius a. On the contrary, in the second mode, denoted as the "static mode", the 
reaction takes place with a finite rate fc, but the sear cher i s imm obile during search phases . Note that this descrip tion 
is commonly adopted in reaction-diffusion systems ( Ricd . [l985l ) or operational research ( Frost and Stond . l200ll ). A 
more realistic description is obtained by combining both modes and considering a diffusive searcher with diffusion 
coefficient D and finite reaction rate k. In order to reduce the number of parameters and to extract the main features 
of each mode, we study them separately by taking successively the limits fc — ^ oo and D oi this general case. 
More precisely, in these two limiting cases, we address the following questions : what is the mean time it takes the 
searcher to find a target? Can this search time be minimized ? And if so for which values of the average durations Ti 
of each phase? 



3. Basic equations 

We now present the basic equations combining the two search modes introduced above in the case of a point- 
like target centered in a spherical domain of radius b with reflexive boundary. Note that this geometry mimics both 
relevant situations of a single target and of infinitely many regularly spaced no n revisitable ta rgets. As in the previous 
model, the mean first passage time to a target satisfies backward equations ( Redneil . l200l[ ) (see section IIV.B.2p for 
derivation), which write here : 

72. 



Dyth + - / - h)de^ - klair)h ^ -1 (14) 

V-Vrt2--(t2-tl) = -l (15) 

where ti stands for the mean first-passage time starting from state 1 at position r, and t2 for the mean first passage 
time starting from state 2 at position r with velocity V, of direction characterized by the angle 6'v- Here Ia(i") = 1 if 
|r| < a and Ia(i") = if |r| > a. In the present form, these integro-differential equations do not seem to allow for an 
exact resolution with standard methods. We thus resort to an approximate decoupling scheme, which relies on the 
following idea. If the searcher initially starts in phase 2, and if the target is close, its initial direction matters. But as 
soon as the initial position is far from the target, there are numerous reorientations before finding the target, implying 
that the initial direction does not matter. Consequently, if & 3> a and once the mean search time has been averaged 
over the starting position, the effect of the initial direction can be neglected. It allows us to ma ke an approximation 
and solve the system (for more technical details, see appendix and iBenichou et aL j (|2006bl [20071) '). 



4. Results for the diffusive mode of detection 

For the diffusive mode of detection {k — > oo), an analytical approximation for the search time can be obtained 
(see appendix IA.2.bp . In the case of low target density (a ^ 6), which is the most relevant for hidden target 
search problems, three regimes arise. In the first regime a <^ b <^ D/V, the relocating phases are not efficient and 
intermittence is useless. In the second regime a ^ D/V <C b, it can be shown (see section rA.2.bp that the intermittence 
can significantly speed up the search (typically by a factor 2), but that it does not change the order of magnitude of 
the search time. On the contrary, in the last regime D/V <^ a <^ b, the optimal strategy, obtained for 

, ^ D \n\b/a) , ^ ^(ln(fe/a) - 1/2) (16) 

^ 21/2 21n(Va) - 1^ y\ \ / J / J ^ V 7 

leads to a search time arbitrarily smaller than the non intermittent search time when V ^ oo. Note that this optimal 
strategy corresponds to a scaling law 



opt 



D 1 



opt n'2 /r, 1 / 1/1 



i2-l/Hb/a)r 



(17) 
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which does not depend on V. 

5. Results for the static mode of detection 

We now turn to the static mode {D — > 0) (see appendix IA.2.a|) . In this case, intermittence is trivially necessary to 
find the target, and the optimization of the search time leads for 6 3> a to : 

/ a X 1/2/ 2 ln(6/a)-i y/^ 

r2,min = ^ (Hb/a) - 1/2)1/' ^ ^^g^ 
which corresponds to the scaling law T2^min = '^krf which still does not depend on V. 

6. Conclusion 

This bidimensional two-state model of search processes for non-revisitable targets closely relies on the experimen- 
tally observed intermittent strategies adopted by foraging animals. Using a decoupling approximation numerically 
validated, it can be analytically solved, allowing us to draw several conclusions, (i) The mean search time (t) presents 
a global minimum for finite values of the r,, which means that intermittent strategies constitute optimal strategies, 
as opposed to Levy walks which are optimal only for revisitable targets, (ii) The optimal r°^* obtained for two modes 
of detection are different and depend explicitly on D and fc, leading to different scaling laws which are susceptible to 
discriminate between the two search modes . (iii) A very striking and non intuitive feature is that both modes of search 
studied lead to the same optimal value of r2^*. As this optimal time does not depend on the specific characteristics D 
and k of the search mode, it seems to constitute a general property of intermittent search strategies. The robustness 
of these conclusions will be discussed further in section IIVI in the framework of a more general model. 



D. Should foraging animals really adopt Levy strategies ? 

As seen before, intermittent strategies are an alternative to Levy walks (defined in section lILAp for interpreting 
trajectories of foraging animals. However, the Levy walks are often thought as optimal and widespread in nature. Is 
it really true? 



1. The albatross story 

Many foraging animals, including albatrosses, deers and b umblebees to name a few, have been thought for long 
to adopt Levy strategies described in the pioneering work of IViswanathan et al. I (I1999D . These foraging behaviors 



were repeatedly accounted for by a simple model stating in the more general framework of search processes that 
Levy walks are optimal search strategies, as they constitute the best way to explore space. Recently, lEdwards et al.l 
( 20071 ) reanalyzed these data, completed by newly gathered data on foraging albatrosses, and showed that in fact 



there was no experimental evidence for the Levy walk behavior^. This study questions the interpretation of several 
experimental works, but also raises a new important and puzzling question : why animals do not adopt the Levy 
walk strategy which has however been reported to be an optimal search strategy? Here we clarify this apparently 
paradoxical situation. 



^ Albatrosses behavior was followed by a humidity sensor on the birds. Flights were taken as the "dry" phases, interspersed with humid 
phases, when the birds touches the ocean. Very long "flights" eventually proved to be rest time, when the bird was in its nest. Once 
these misinterpreted dry phases removed, the distribution of flights durations is no longer a power-law. 
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2. Do animals really perform Levy walks? 



As the optimality of Levy strategies crucially requires conditions on the targets (regenerating at the same place, 
patched or hard to capture) and conditions on the searcher (no switch when a target is found, which is a very simple 
form of memory), it can not be taken as a general rule even if realistic for certain species. On the contrary, we argue 
that the general question of determining the best strategy for f inding a single hidd e n targ et belongs to the situation 
of destructive search, where in the framework of the model of IViswanathan et al.l ( 19991 ). the most efficient way to 
find a randomly hidden target is simply a linear ballistic motion and not a Levy strategy fs ee section III. All . As a 
consequence, there is no paradox : the reason why Levy walks are not observed in the work of lEdwards et al.l (|2007h 
is probably because they do not constitute robust optimal search strategies. 

And what about other experimental observations? Among experimental studies analyzing organisms trajectories 
as a succession of segments interspersed with turns, an important prop ortion reports t i mes between turns distributed 
expo nentiall y (a list of e xamp les far from exhaustive : C.Elegans worm (iFuiiwara et all. 120021: [Pi^rce-Shimonura ej^L 
Il999l). fish (iHill et al.l . |2000|) . plankton in a part of the conditions studied bv iBartumeus et al.l (I2003D . amoebae 



(iLi et all I200S 



However, apart from the controversial albatrosses study (Edwards et al., 12007 ). there is 



boom in articles claiming that Levy behavior is observed for some anim al species. A pa r t of t hem can be dismissed 
as evidence of Levy behavior. On the one hand, as explained in details in lEdwards et al.l ( 200l1 ). due to experimental 
limitations, most data cover only a very limited range, which makes difficult a reliable identification of power laws. On 
the other hand, patterns and processes should not confused, as underlined bv iBc nhamou (2007i ). The same observed 
patterns, can often be explained by different models. It is not because a trajectory is similar to Levy walk trajectories 
that the underlying process is necessarily a Levy walk. For example, a composite classical random walk can look very 
similar to a Levy walk for a time short enough (see figure [7]). Nonetheless, not all studies should be discarded , since 
li mited studies neither prove nor rule out Levy strategies, fsee IViswanathan et al.l ( 2008f l for a review). As underlined 
bv IViswanathan et al.l ( 20081 ). other selection pressures could be predominant. For example, when targets location is 
known, exploitation could be opti mized instead of search, and Levy walks could emerge fr om interactions between 
the environment and the searcher ( Bover et al.l , l2006l : iJiang et al.l , l2009t ISantos et al.l , l2007t ) . 



E. Conclusion on animal foraging 

Levy walks are a fashionable model for interpreting trajectories of foraging animals. However, on the one hand, 
there is controversy about at least a part of the experimental data which were thought to support Levy walks. On 
the other hand, the conditions in which Levy walks are optimal are very restrictive. However, it do es not rule out 
any co ntribution of Levy statistics in the context of search processes. For example, as discussed by iLomholt et al.l 
(I2008D and in section IV.C.21 Levy statistics can be advantageously used in the context of intermittent trajectories. 
Additionally, we argue that some animals cannot detect their target when they are moving ballistically, and in fact 
alternate these fast but blind phases with detection phases. Mean search time with intermittence can be smaller than 
with detection phase alone, and it can be minimized tuning the mean durations of each phase. Intermittent search 
strategies, because they rely on the experimental observation that speed degrades perception, and because they prove 
optimal an robust, are good candidates for interpreting animals trajectories. 



III. INTERMITTENT SEARCH STRATEGIES AT THE MICROSCOPIC SCALE 

The previous section has shown that intermittent search strategies are observed at the macroscopic scale. They are 
also observed at the microscopic scale. In the following, we will focus on two examples : the localization by a protein 
of a specific DNA sequence, and the active transport of vesicles in cells. 



A. Protein/DNA interactions 

1. Biological context 

Various functions of living cells - and therefore at larger scales of living organisms - are regulated by coordinated 
chemical reactions between specific molecules, which are often present only in a few copy number. The importance 
of the kinetics of suc h search processes between reaction partners can be illustrated by the bacterial restriction and 
modification system ( Wilson and Murravl [T99ll ) . which involves couples of methyltransferase an d restriction enzymes 



that recognize the same sequence on DNA (for example iJcoRV recognizes the sequence GATATC (|Tavlor and Halfordl . 



1989D ). Methyltransferase enzymes methylate this specific sequence on the bacterial DNA in order to protect it from 
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Levy walk with jx = 2 



Composite random walk : alternation of about 10 short steps (mean 1, distributed exponentially), and one large step (mean 

10, distributed exponentially). 

Figure 7 Comparison between a Levy walk and a composite random walk : They are not easy to distinguish at short time 
scales 



restriction enzymes, whose function is oppositely to cut the DNA at this specific sequence. This function is first aimed 
at impairing any intruder viral DNA which enters the cell and which is very likely to contain the target sequence. 
Indeed, this sequence, typically 4-8 base pairs, is very short as compared to the viral genome, which, depending on 
the virus, can be made of lO'^ — 10® base pairs (typically 5.10^ for bacteriophages). The infected bacteria then faces 
a vital search problem : restriction enzymes must find their target sequence on the viral DNA reliably to inactivate 
the virus before it exploits the bacteria machinery and kills it. 

More generally, it is well established that so me sequence - specif ic proteins find their target site in a remarkable 
short time. For the lac repressor for example, iRiggs et al.l ( 1970l ) measured association rates orders of magnitude 
larger than those expected for reactions limited by the classical three-dimensional diffusion (results confirmed by 
Hsieh and Brenowitj ( ToOVt) at different salt concentrations, ruling out electrostatic effects as the only explanation). 
HalfordI (|2009ir argues that in fact only a few enzymes react significantly faster than the 3-dimensional diffusion limit. 



However, this author underlines that many enzymes react at rates close to the diffusion limit, and that this observation 
is still impressive. Indeed, classical experiments are performed with a considerable excess of DNA, which is very likely 
to contain sequences similar to the target sequen ce and therefore act as t raps slowing down the enzymes in their search. 
In a series of seminal articles, iBerg et all ( 198l[ ); IWinter et al.l ( 198l[ ): IWinter and Von Hippell ( 198l[ ) proposed that 
3D diffusion (or "hopping"/"jumping") was not the only motion available to the protein, even if no energy is consumed 
(unlike some enzymes which consume energy to scan DNA processively) . They suggested that in some cases, proteins 
could bind non sp ecifically to DNA due to a weak elec t rostat ic interaction and diffuse along the chain in a process 
named sliding (see I Von Hippell ( 20071) and lDahirel et al" jilooi) for more details on the weak electrostatic interaction). 
It was then argued that the combination of sliding and 3D diffusion, i.e. facilitated diffusion, can make the search 
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Figure 8 Artist view of a DNA/protein interaction, which combines 1-dimensionaI sliding phases and 3-dimensional relocation 
phases. Picture by Virginie Denis/Pour la Science 352, February 2007. 



for a sequence two orders of magnitude faster than 3D diffusion alone and henceforth sufficiently efficient (see also 
lAdam and Delbruckl (|l968l )). 

Actually, this search mechanism can be classified as intermittent, in the general meaning defined in introduction. 
Indeed, on the one hand, 3-dimensional diffusion off the DNA is fast, but does not allow for target detection. On 
the other hand, sliding is a phase of motion along DNA, which therefore enables target detection, but which is much 
slower due to a higher effective friction. 

The pioneering studies on facilitated diffusion ( Berg et al.l . llQSlt iRiggs et al. I Il970l : IWinter et al. 
IWinter and Von Hippell . Il98lh are based on ensemble measurements, which were for a long time the only way to 
experimentally access to protein/DNA interaction. Recently developed techniques make possible the observation of 
this interaction at the level of a single molecule, with a resolution in space a.nd tim e still improving (for a recent 
review on the experimental results of such techniques, see iGorman and Greend (120081) '). ft is now confirmed directly 
that many proteins searching for a specific sequence on D NA combine "hoppin g /jumping" a nd "sliding" (see figure [HI) . 



Sliding phases have been clearly id entified (both in vitro ( Kabata et al.l . Il993[) and in viv o ( Bak. 



Elf et al.j. I2OO7'; 'Wang et_al.', '2006')) , as well as hopping /j umping phases ( Bonnet et al.l , l2008t 



and Metzler 



Gowers et al 



2004 : 



2005; 



Komazin-Meredith et al.. . ,2008; .van den Broek et al. [200g 

With these new single-molecule experiments, theoretical models have bloomed too. Here, we first present a 
stochastic approach of a simplified version of the problem, which shows that it is the intermittent nature of the 
trajectories that makes possible such high reaction rates. This minimal model permits to calculate explicitly the 
mean search time for such intermittent reaction paths, and shows that reactivity can even be optimized by properly 
tuning simple dynamic parameters of intermittent trajectories. Next, we discuss the different directions of extension 
of recent theoretical models. . 



2. Minimal model of intermittent reactions paths 



We present here a simple model of intermittent reaction paths with minimal ingredients ( Coppev et al.l . [2004[ l. We 



first define the model, then explain the main steps of the calculation, and eventually give the results. 
Definition of the model. 

We consider a generic protein searching for its target site on a DNA molecule (see figure |9]) . The pathway followed 
by the protein, considered as a point-like particle, is a succession of ID diffusions along the DNA strand (sliding 
phases denoted phases 1) and 3D excursions in the surrounding solution (denoted phases 2, during which the target is 



17 



protein 




DNA— ^ 



Figure 9 A model of intermittent transport for DN A/protein interactions. 

not accessible). In this minimal model we assume that the target site is a perfect reactive point of the DNA molecule, 
which means that reaction occurs as soon as the target is reached by the protein. Note that in this case the protein 
can find the target site only by diffusing along DNA, and therefore follows the scheme of intermittent search presented 
in introduction. The ke y quantity that we evaluate in this section is the search time, or reactio n time, defined as the 
mean first-passage time ([Benichou and Voituriej . l2008t[Condamin et al. I l2007al l2005allbl . I2007hll3 . l2008t iRedne j . 1200 iD 
of the protein at the target, here denoted by ( T). This quantity gives dir ectly access in a mean field approximation 
to the first order reaction constant K ~ 1/ (T) (|Berg and Blomberel Il976l ) . 

We now introduce further ingredients of the model. The time spent by the protein on DNA during each sliding phase 
is assumed to follow an exponential law with dissociation frequency Ai. This law relies on a Markovian description 
of the chemical bond which is commonly used for bimolecular association. The probability density that the protein 
leaves the DNA at a random time t is then given by Ai exp(— Ait)dt, and the mean duration of a sliding event reads 

Tl = 1/Ai. 

The one-dimensional motion on DNA during sliding phases 1 is modeled by a continuous Brownian motion with 
diffusion co efficient D. We assume that the ends of the DNA chain act on the protein as reflecting boundaries (see 
for instance iJeltsch and PingoudI (|l998f )). but in practice this assumption is unimportant for long DNA molecules. 



Moreover, the case of circular DNA, such as plasmids, is readily obtained by taking the particular case L = R. We 
next assume that the 3D excursions of phase 2 are uncorrelated in space. This means that after dissociation from 
DNA, the protein will rebind the DNA at a random position independently of its starting position. This is justified 
when the DNA is in a random coil conformation, as in this case even short 3D excursions can lead to a long effective 
translocation of the linear position of the protein on DNA. We further assume that the probability density Paoit) of 
the duration t of such 3D excursions is exponentially distributed, and write Psoit) = X2 cxp(— A2t). This assumption 
is justified, at least for the tail of the distribution, as soon as the 3D excursions of the protein are confined in a closed 
volume, for instance an experimental volume in vitro, or in vivo the cell or a cell compartment. The mean time 
T2 = I/A2 spent in the surrounding solution in phase 2 can then be shown to be proportional to the c onfining volume 
(jBenichou et"aH , l2005al . l2008t iBlanco and Fournieil . 120031 ICondamin eFaD , l2005aL l2007bl : [k3, [Hli). 



We next introduce PiD{t\x), which is the conditional probability density that the protein, being on the DNA at 
position X and at time t ~ 0, will dissociate at time t before any encounter with the target site. We rewrite this 
quantity as : 

PiDit\x) ^ Al cxp(-Ait)Q(t|a-) (20) 

where Q(t\x) is the conditional probability density that the protein, starting from the position x, does not meet the 
target site during a single sliding event. The probability density j(i|a;) of the first passage to the target site position 
at time t without dissociation is then related to Q(t\x^ according to Q(t\x^ = 1 — j(^'|a;)dt'. 

Last, we introduce Pi£)(<|a;), which is the conditional probability density that the protein, being on DNA at position 
X at time t = 0, will find the target site for the first time at time t within a single sliding phase 1, without leaving 
the DNA: 

Pii3(ik) =exp(-Ait).7(iN). (21) 



Given these quantities, we show below that the first-passage density of the protein to the target site, and conse- 
quently the reaction constant, can be calculated explicitly. 
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First-passage density. 

By calculating the first-passage density, we obtain the mean reaction time, as well as all associated moments. We 
assume that the protein starts at t = in state 1 (bound to the DNA) at position x. We consider a generic event 
whose number of 3D excursions is ti — 1, and denote the duration of successive sliding phases ti, . . . ,t„, and the 
duration of successive 3D excursions 6i,. .. ,On-i- The probability density of such an event, for which the protein 
finds the target site for the first time at time t = X^ILi + Y^7=i ^'^^ 



Pn{t\x) = PlD{tn)PMSn~l)PlD{tn) ■ ■ ■PlD{t2)PMSl)PlD{tl\x) 



(22) 



where Pioit) and Pmit) are averaged over the initial position of the protein: Pioit) = {PiD{t\x)) ^ and Pmit) — 
{PiD{t\x)') ^. We denote by L the DNA length on the "left"' side of the target site and by R the length on the "right" 

side of the target site. The average of a function / over the initial position x is given by {fit\x))^ = J^j^ f{t\x)dx. 

To obtain the density of first passage to the target site, F{t\x), we sum over all possible numbers of excursions and 
we integrate over all intervals of time, ensuring that t = ti + 9i. The average over the initial position of the 

protein, F{t) = {F{t\x))^, can be expressed as: 



F{t) = XI / • • • '^*nd0i . . . dOn-lS [^U + ^d,-t 



(23) 



n P^D{e^) 



1=1 



n p^d{u] 



Pioitn). 



Taking the Laplace transform of F{t), F{s) = dte ^*F{t), we obtain: 



Fis) = (j(X,+s\x) 



j(Ai + s\x] 



(l + s/Ai)(l + sA2) 



(24) 



where j{s\x) is the Laplace transform of j(t\x). This expression completely solves the problem for any ID motion. 
We will see next that the main quantities of physical interest can be extracted from this formula. 

Optimal search strategy. 

The relevant quantity to describe the protein/DNA association reaction is the mean time (T) necessary for the 
protein to find the target site (see above) . This mean time is obtained from the derivative of the first passage density 
by the following relation: 



which combined with equation (j24p gives: 



(T) 



dFjs) 
ds 



(25) 



s=0 



(T) 



l-(j(Ai|x))^ 



.?(Ai|a:; 



1 



1 



(26) 



This expression is general and holds for any ID motion in the slow phase 1. Now, we calculate this quantity in the 
case where the phase 1 is a free ID diffu sion. The one-dimensional Laplace transform of the first passage probability 
density is well known (see iRediieil ( 200ll) ). and leads to, after averaging over the starting position x: 



^{L + R) 



tanh ( J 



tanh ( \/^R 



- 1 



(27) 
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where D is the diffusion coefficient. This defines as a by product the association constant of the reaction as K = 1/ (T) . 
Two first comments are in order, (i) First, as soon as the length of the DNA strand is large enough (more precisely 

as soon as \J^L ^ 1 or ^J^R ^ 1). (T) grows linearly with the length of the DNA. This mirrors the efficiency of 

intermittent reaction paths, as compared to the quadratic growth obtained in the case of pure sliding. In particular, 
the boundary effects are negligible for this quantity as soon as the overall length is large enough, (ii) Second, this 
expression is valid for a very large class of 3D motions. More precisely, it holds as soon as the mean first return time 
T2 corresponding to the 3D motion is finite and independent of the departure and a rrival points. 

We now come to an important question recently addressed in lCoppev et al 1 (|2004 ): ISlu"tskv and Mir^ (|2004h . which 
concerns the optimization of such intermittent reaction paths. We here assume that the mean search time is a limiting 
quantity which might have been minimized in the course of evolution. In this context, we consider Ai = 1/ti, which 
characterizes the protein/DNA affinity, as the adjustable parameter. Indeed, this quantity depends strongly both on 
the structure of the protein and on physiological conditions such as the ionic strength, and therefore could widely 
vary from one protein to another. On the contrary A2 = I/T2 depends mostly on the properties of the environment, 
such as the DNA conformation, which is itself subject to very stringent constraints and therefore much less likely to 
be varied. Another adjustable parameter is the ID diffusion coefficient D. Optimizing the search time with respect to 
this parameter is trivial: it is found that D should be as large as possible (assuming that D and Ai are independent), 
but obviously one should keep in mind that D is controlled by the hydrodynamic radius of the protein which can not 
be too small. For these reasons we focus here on Ai. 



(T) 



diff 
5- 



0.5 



0.1 




10" 



10" 



10" 10^ 

n/T2 



10^ 



10" 



Figure 10 Mean first passage time of a DNA-binding protein to its target site, renormalized by the mean first passage time 
by diffusion alone (Tdiff), as a function of ti/t2. Tdij] Iti — 1 (red dashed-dotted line), Tdiff/T2 = 10 (violet dashed line), 
Tdiff/T2 = 100 (blue dotted line), Tdiff /t2 ~ 1000 (green solid line). For a small DNA length £ (and therefore small Tdiff) 
(red dashed-dotted line), the reaction constant depends monotonously on ri and intermittent reaction paths are inefficient. 
For larger values of I (other curves), the reaction rate can be optimized as a function of ri. Here we averaged other the initial 
position of the target. 



It can be seen qualitatively that (T) is large for both Ai very large (in the limit Ai infinite, the protein is never 
on the DNA), and Ai very small (pure sliding limit which gives a quadratic growth with the DNA length), and could 
therefore be minimized for and intermediate value of Ai. The sign of the derivative of the mean search time at Ai = 
shows that it can indeed be minimized provided that 

A2 > 15D , (28) 

This condition means that bulk excursions, to be favorable, should be shorter than a fraction of the typical time 
needed to scan the full DNA molecule by ID diffusion. In particular, it requires that the DNA length is long enough. 
If this condition is fulfilled, the search time can indeed be minimized (see figure [TU)) . A careful analysis of the implicit 
equation satisfied by Ai at the optimum leads to the following expansion for large £ = L + R 
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Equations ((28|) and p9)) refine the result of ISlutskv and Mirnvl ( 2004 ). which predicts than the optimal strateg; 



is realized when Ai = A2. This result actually holds in the large i limit, or more precisely for \J ^£ ^ 1. For 
intermediate values of £ boundary effects become important and the minimum can be significantly different. 

The (T) value at the minimum is particularly interesting. We compare it to the case of pure sliding where (T,) ~ 
£^/i3D): 

The efficiency of the 3D mediated strategy is therefore much more important when the DNA chain is long. For 
example, using standard values for Ai (a few 10~^ s.) and D (typically 10'^ A*^/s-) and for a DNA substrate of length 
10^ bp, the mean reaction time is three orders of magnitude smaller than for a pure sliding strategy. Beyond the 
importance of such results for understanding the kinetics of gene transcription, this first minimal model shows that 
intermittent reactive paths are indeed very efficient, and that they can even allow to optimize the reaction kinetics. 



3. Towards a more realistic modeling 

The model introduced above provides a simple way to discuss the minimization of the search time. Further 
approaches have been developed to model target search by proteins. We present below the main models used in 
the literature, and discuss their relevance to real target search problems by proteins in cellular conditions. 



a. Main approaches 

Generally speaking, theoretical models of facilitated diffusion rely on the basic assumption that the protein alternates 
phases of ID diffusion along the DNA and phases of free diffusion when the protein is desorbed from the DNA. The 
existence of such two distinct states, whose dynamics is usually characterized by association/dissociation rates, is 
supported by direct experimental observations as discussed above. Additionally, molecular dynamics simulations 
taking into account the electro static interact i on be t ween the negatively c harge d DNA and the locally positively 
charged protein (see for example iDahirel et al.l ( 20091 ) : lFlorescu and Joveip] ( 20091 )) have shown that these two states 



naturally arise on the basis of the electrostatic interaction only, suggesting the robustness of the facilitated diffusion 
mechanism. Such studies at the molecular scale could serve as a tool to calculate the association/dissociation rates 
used in the models of facilitated diffusion discussed in this section, which all take into account only effectively the 
molecular interactions. 

Stochastic modeling 

The minimal model presented in section lIII.A.2l relies on the statistical analysis of the trajectory of a single protein, 
and can henceforth be qualified as a stochastic model. Similar stochastic methods have been used and complemented 
inlBinichou et al. (2009); Efiazar et al, (2001., .2008, ): .Lomholt et al. (2005, 2009, 2007): Mero^ ct_al. (2003), ■•■ and 
have the advantage, when solvable, of giving access to the full distribution of the search time, yielding a refined 
information on the search kinetics. Moreover they can be adapted in some cases to take into account anomalous 
transport both in the ID and 3D phases as discussed below. 

Kinetic approach 

The main alternative to the stochastic approach is given by what can be called kinetic models, which assume a 
steady-state homogeneous concentration of proteins in contrast with the single protein description of stochastic models. 
Such models therefore rely on a mean-field approximation which proves to be efficient to evaluate the mean search 
time thanks to scaling arguments. A first example is given bv lHalford and Markol (|2004[ ). where scaling arguments are 



used to roughly estimate the time for the protein to find the DNA coil, and then the time to find the target inside the 
coil, whic h eventually yields an optimal sl i ding l e ngth. More generally, the key ingredient of kinetic models, developed 
mainly in lHu L et al. I (|2008^ : lHu T et all (I2006D :[h u and Shklovskiil (j2006 ) , is that the system is assumed to be in a 



stationary state. Under this hypothesis, the fiux of particles delivered by the 3D diffusion into the sphere of influence 
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of the target, whose size is defined as the "antenna length" ^a, must be equal to the flux of particles delivered by ID 
diffusion into the target. Such balance equation generically reads 

J ^ D'iCfreeia ~ DlCads/K, (31) 

where the concentrations of free (c/ree) and adsorbed (cads) proteins are assumed to be at equilibrium, ie satisfying 
Cf reel Cads = K whcrc K is the equilibrium constant associated to the association/dissociation rates. Importantly, the 
antenna length, defined as the typical scale below which the dominant transport is sliding instead of 3D diffusion, has 
size when measured in 3D space, but takes another value \a when measured along the DNA. Making assumptions 
on the DNA conformation (for instance random coil or fractal globule), different scaling laws between Aa and can 
be proposed. Equation ([5T|) then permits to determine and henceforth to give the scaling of the mean search time 
1/J. The advantage of this method is that it permits, through the relation between A^ and ^q, to take into account 
various models of DNA conformations, which is much harder to achieve in the stochastic approach. Such models, 
whose results are compatible with the stochastic approach, provide in addition a useful picture of facilitated diffusion. 
Indeed, in these models the effect of sliding can be seen as effectively making the target of the size of the antenna 
length, which is much larger than the real target size, and therefore speeding up the search. 



Finally, these two approaches are quite complementary and both require as an input the modelling of ID and 3D 
phases. The minimal model of section ITlI . A. 2 1 describes the ID phase as regular diffusion, while 3D phases are assumed 
to result in completely random relocations over the DNA. Beyond this minimal model, the specific description of these 
two phases has motivated numerous works and many refinements have been discussed in the literature. We review 
in the next sections the main models that have been proposed to provide a more realistic description of ID and 3D 
phases. 



b. Descriptions of the ID phase (sliding and recognition ) 



diffusion, 
or otiier 
motions 




Figure 11 Sliding is often represented by diffusion with perfect reactivity on the target. This figure shows the two main 
directions for a more realistic description of the ID phase : on the one hand, the shding is not necessarily diflFusive, and on the 
other hand the ID phase could be in fact a combination of 2 phases, one fast but with low recognition, and another slow (or 
immobile), but with high recognition of the target. 



Anomalous diffusion in the sliding phase 

As stated above, the phase of one-dimensional non specific interaction of a protein with DNA, "sliding", is generally 
described as Brownian di ffusion as in the minimal model of sectio n IIII.A.2I If this hypothesis seems to be confirmed 
by in vitro experiments ( Bonnet et al.l , l2008l : iKabata et al.l , Il993l ). it could not always be the case, in particular in 
vivo. A first limitation of this simple description appears in the case of many protein s binding to DNA, as it is the 
case in vivo, which are likely to create traffic jams ( Li et al.l . 120091 : ISokolov et al.Ll2005l ). Such crowding effects in one 
dimension are known to potentially lead to a subdiffusive behavior. Additionally, even in the case of a single protein, it 
should be kept in mind that the DNA sequence is not homogeneous, and the disorder in the sequence can also impact 
on sliding. The heterog eneity i n the sequence is o ften modelle d by a disordered energy landscape, whose distribution 
is Gaussian (jBarbi et al. , 2 004. : ,Hu and Shklovskiit . .2006;.Wunderlich and Mirnvl . l2008h . iBarbi et al. (20041 show that 
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in this case sliding is not purely diffusive : at short times, the protein will be trapped in local minima, leading to 
a subdiffusive behavior. The diffusive behavior is recovered only at larger times, or equivalently for sliding lengths 
longer than a hundred of base pairs. 

Anom alous diffusion i n the sliding phase has been discussed at the theoretical level in lEliazar et al. I (I2007D fsee 



also.Eliazar et al. 1 (l2008^ : [MCTOz et all (I2009D ). which extend the minimal model of facilitated diffusion ICoppev et al.l 



( 2004 ) summarized in section IIII.A.2I In particular, the Laplace transformed search time distribution is obtained 
for several non Brownian sliding motions such as ballistic, self-similar or halted motions (in particular when halts 
durations are widely di stributed, leadhrg to a subdiffusive behavior) . therefore covering standard models of anomalous 



diffusion. Importantly, lEliazar et al. find that whatever the model of sliding, there are always regimes in which 



intermittence is favorable, similarly to the case of Brownian sliding. They further show that in the case of 3D 
excursions with finite mean durations, the mean search time with an arbitrary sliding mechanism remains of order 
oc £. This indicates that for long enough DNA, intermittence is favorable for a wide range of sliding motions, be it 
normal or anomalous, which supports the robustness of the facilitated diffusion mechanism. 

Target recognition 

A simplification which is often used in the literature consists in assuming that the target is perfectly reactive, 
i.e. that reaction occur s with probability one at the very first passage of the protein to the target sequence. 
ISlutskv and Mirnvl (|2004t ) however stress that if there is an activation barrier at the target site, the protein as a 



chance to pass on the target without entering the recognition process, and therefore to m iss it. The roughne s s a o f 
the energy landscape of the sequence can then play a crucial role. If a is of order IfcsT, ISlutskv and Mirnvl (|2004l ) 



find that the sliding diffusion is fast, but that the protein has a high chance to miss the target. On the contrary, if a 
is of order bk^T, the recognition probability is high, but the sliding diffusion coefficient is very low, leading to a huge 
search time. Such result seems to set confiicting constraints on the search, sinc e an efficient target search process 
requires both speed and reliability in target rec ognition. To overcome t his paradox. ISlutskv and Mirnvl ( 20041 ) propose 



on the basis of direct structural observations ( Kalodimos et al.l l2004t ) that the protein can perform conformational 



changes in the sliding phase, and switches between a fast search state (with low a) and a slow recognition state (with 
high a). They show that if the two energy landscapes are strongly correlated, it is possible to conciliate high speed 
and efficient recognition. Similar models of 2-state proteins in the sliding phase have been stu died more quant i tative ly 
both in the framework of the kinetic (see lHu L et al.l (200l)) and the stochastic approach (see lBenichou et all ( 20091 )). 



confirming that such mechanism indeed permits a fast search with reliable recognition. 



c. Descriptions of the 3D phase (jumping/hopping) 

We now review the different descriptions of the 3D phase. Most of models of facilitated diffusion implicitly require 
the knowledge of the probability 7r(x|a;o) that a protein which desorbs the DNA at position xq will eventually rebind for 
the first time to the DNA at position x, where the coordinate x measures the distance along the chain. This quantity 
depends both on the dynamics of the protein and on the conformation of the DNA, and in practice can be determined 
explicitly only in the case of an ideal infinite cylindrical DNA (see below), which makes assumptions necessary in 
realistic situations. Depending on the relocation length \x — xo\, 3D excursions have been given in the literature 
different names, mainly either "hops" (referring to "small" \x — xo\) and jumps (referring to "large" \x — xo\). Since 
the definition of jumps and hops may vary according to authors, the limit between both being somehow arbitrary, we 
give below the one that will be used in this review. 



Jumps 

We define as jumps the 3D excursions whose starting and ending points on the DNA sequence are uncorrelated, 
ie such that the relocation probability 7r(a;|xo) is independent of x. This definition of course depends on the DNA 
conformation, and makes sense only for dense enough conformations such as a random coil (for insta nce for free DNA 
in so lution) or even denser packings that can be expected in vivo (such as fractal globule structures ( Grosberg et al.L 
Il988l ) ). In such structures, sequences which are far apart along the DNA chain can be actually very close in the 3D 
space. Hence, 3D excursions whose relocation length is larger than the typical distance between DNA segments are 
likely to end at any remote location on the DNA sequence, and should therefore be considered as jumps according to 
our definition. The lower bound of the relocation length of jumps therefore strongly depends on the DNA conformation 
and is in practice hard to evaluate. In the case of interest of a confined DNA (as in cellular conditions), the typical 
distance between DNA segments can be estimated by the DNA density-density correlation length ^c- This suggests 
an alternative and equivalent definition of jumps as 3D excursions whose relocation length is larger than ^r.. A 
widespread assumption used for exam ple in the minimal mo d el of section IIII.A.2I ( see also Hi^Te^L (|2006l ) at 
scales larger then the antenna size or iBenichou et aL I (|2009D : ICoppev et al.. (,2004! ): lEhazar et all (j2007r )) consists 
in taking all 3D excursions as jumps, ie as random uniform relocations over the DNA. Even if not exact, this 
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Figure 12 Description of 3D excursions. Jumps are 3D excursions whose starting and ending points on tlie DNA sequence are 
uncorrelated. In practice, for confined DNA conformations as in cellular medium, jumps have a span (measured along the DNA 
contour) larger than the density-density correlation length of DNA. Conversely, hops are 3D excursions whose starting and 
ending points on the DNA sequence are correlated, or equivalently 3D excursions whose relocation length is smaller than the 
correlation length ^c- 



assumption has been checked numerical ly on the example of a quenched self avoiding DNA and as proved to be very 
satisfactory ( Sheinman and Kafril 120091 ) at high enough DNA concentration. Interestingly, this defi nition of jumps 



highli g hts the importance o f the l ocal DNA concentration on the search efRciency, as observed in iLomholt at al.l 
(|2009D :I van den Broek et al.l (j2008l ) : the more densely packed the DNA, the smaller the correlation length ^c, and 



therefore the higher the probability of jumps enabling the protein to explore previously unscanned areas, and the less 
the time spent in 3D phases. 

The assumption that all excursions of relocation length larger than the typical distance between DNA segments 
lead to a completely uniform relocation over the DNA is however not exact. Indeed, large scale correlations in the 
3D conformation of the DNA chain may exist, as in the model case of a free random coil conformatiq r i. In p articular 
in eukaryots, where DNA is packed in the nucleus, recent studies such as ( Lieberman-Aiden et al.l, |2009| ) support 



a hierarchized structure of DNA which could induce long range correlations in the conformation and therefore a 
non uniform relocation probability. The impact of DNA con formation on the relocation p robability has been tested 
on the example of an annealed worm-like chain polymer by iDiaz de la Rosa et al.l ( 2010f ) who found that at short 



times, 3-dimensional excursions of approximately the DNA persistence length are actually less abundant than both 
shorter and a bit larger relocations : by definition, closer sequences along the DNA are closer in 3D space, but 
sequences further than the persistence length can take advantage of loops and actually be even closer in 3D space. 
More generally, the l oop statistics impacts on the relocation probability and can lead to a non uniform relocation 



probability 7r(a;|a;o ) . ILomholt et all (|2005l ) argue that since polymers form loops whose linear size x is distributed 



according to p{x) ~ |a;|~ (for instance a = 0.5 for Gaussian chains, a ~ 1.2 for self-avoiding walks), relocations 
distributed according to the same law are favored and should also be taken into account. In an annealed version of 
such a model, the authors show that depending on a the optimal strategy can widely vary. It should be added that 
loop s can enable another reloc ation mechanism for prote i ns wit h multiple binding sites called intersegmental transfer 
(seeE u and Shklovskiil (|2003) and |Sh einman and Kafril (|2009[ )) which can be shown for modelling purposes to be 



widely equivalent to 3D excursions. Finally, the approximation of uniform random relocat ions proved to be useful 
and can be validated numerically for simple DNA conformations ([Sheinman and Kafril2009l) . but a better knowledge 
of the in vivo conformation of DNA would be necessary to assess more precisely the relocation probability. 

Hops 

Echoing the definition of jumps, we define hops as 3D excursions whose starting and ending points on the DNA 
sequence are correlated, or equivalently as 3D excursions whose relocation length is smaller than the typical distance 
between DNA segments, which is given by the correlation length for con fined DNA. Due to their local character, 
hops are often effectively taken into account in the sliding mechanism (see iHu T et all (|2006f) for example). If this 



assumption is in practice useful, it also raises additional questions, (i) First, hops do not explore continuously the 
DNA, and a protein performing hops have higher chances to miss the target (note that hops permit as a counter part 
to bypass obstacles on the DNA), and this non perfect reactivity has to be taken into account, (ii) Second, if hops are 
included in an effective sliding mechanism, then the effective diffusion coefficient as to be determined, as well as the 
effective transition rate from this effective sliding state to the jumping state. This last point amounts in practice in 
calculating the relocation probability 7r(a;|a;o) of hops, which gives as a byproduct the probability that after desorption 
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Figure 13 Facilitated diflFusion of a protein on DNA. Left: schematic definition of sliding, hopping and jumping. Hops are 
3D excursions whose starting and ending points on the DNA sequence are correlated, or equivalently 3D excursions whose 
relocation length is smaller than the typical distance between DNA segments. In the case of confined DNA, this distance is 
estimated by the DNA density-density correlation length ^c. Right: model parameters. 



from DNA, the protein performs a jump rather than a hop. 

This p roblem, which has bee n studied numerically in IWunderlich and Mirnvl ( 2008[ ). can actually be studied ana- 
lytically ( Loverdo et al.l , [2009cl) . since at the scale of hops, which are smaller than the DNA persistence length ^p, 
the DNA can be simply modelled as a cylinder of effective radius a = Rdna + Rprotein (see figure [13] left). Denoting 
D the 3D diffusion coefficient of the protein (assumed point-like). The density probability 7r(a::|ro) of being adsorbed 
on the DNA at the longitudinal ab scissa x, starting frorn the point rp then satisfie s (for d etails on derivation of thi s 
equation and the fo llowing ones, see lLoverdo et all |2009cl ). and lBerg and Blomberd ( 19761 ) : [Chechkin et al. 
iLevitz et al.l (200l) for similar analysis): 



Aro7r(x|ro) = 



(32) 



Assuming radiative boundary condition on the cylinder surface to account for the adsorption rate on DNA, 



droTT{x\ro){ro 



a) = K7r(a;|ro)(ro = a). 



(33) 



one can show that: 
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where Kj ar e Bess el functions. This relocation distribution, in good agreement with experimental data from 
[Bonnet ct al.l obtained on an extended DNA, enables to address the questions of point (ii) above. First, 

Eq. (|34p gives the analytical distribution of hops, since for x < £,c the DNA can be well approximated by a cylinder. In 
turn, as shown in ([Loverdo et al. j. l2009cD . this gives access to the effective diffusion coefficient of a combine motion of 
sliding and hops. Second, according to our definition, all relocations with x > will be jumps. Thus, the probability 
that a 3D relocation is a jump rather that a hop is given by the complementary cumulative distribution: 



C{x = Cc) 



Tr{x\ro)dx 



kl>Cc 



ln(ro/a) + l/na 



(35) 



Going back to the search probl em, in regimes where j umps are favorable, decreasing the correlation length speeds up 
the search process, as found bv lvan den Broek et al.l ( 20081 ). 



Crowding effects 

In both prokaryots and prokaryots, the cellular medium is very crowded, and 3D excursions of proteins are likely 
to be hindered. The normal or anomalous nature of transport in cellular medium is still debated. For instance. 
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iDix and Verkmanl ( 2009( ) supports that 3D motion is mostly normal diffusion; whereas iMalchus and Weisi ( 2010l ) 
suggest that there is more and more evidence for subdiffusion. Experimentally, it is found that some tracers exhibit 
transient behavior (subdiffusive at small timescales, diffusive at larger time scales) : measures however depend on 
the size and nature of the tracer, on the time and length-scales covered and other experimental conditions, which 
may explain the lack of consensus on the problem. The influence of subdiffusion on target search will depend on 
the microscopic mechanism at play. There are three main mechanisms leading to subdiffusion, as underlined section 
[K\ : random walk on a fractal medium, random walk with long waiting times (CTRW), or random walk with long 
range correlations such as Fractional Brownian Motion (FBM). Which of these pos sibilities best de sc ribes transport 
i n crowded environments such as th e cellular medium is ho wever still unclear fsee iHe et al.l (|2008l ). iBancaud et al.l 
(l2009D.ISzvman ski and WeissI (|2009t ) or lTeiedor erall (|2010t ) for various opinions on the subject). 



Lomholt et al. (2007) explore the effect of a crowded environment with subdiffusion {r (t)) oc P (0 < a < 1) caused 



by waiting times distributed as p{t) ~ T"/t^+". They argue that because of these waiting times, the probability that 
the protein has not yet left the DNA at time t and the probability that an unbound protein has not yet bound to 
DNA after a time t both scale as f/t^+". Their results have two main practical implications. On the one hand, in an 
experiment, since proteins can remain stuck for v ery lo ng times, ensemble averages do not lead to the same results 
as time averages as also highlighted bv iHe et al. On the other hand, since proteins would slide for a longer 

time and as it would take them a very long time to ret urn to DNA, the genes coding for transcription factors should 
be close to their target sequences, as also underlined bv IWunderlich and Mirny (12008). The anal ytical determination 
of the relocation distribution above can be extended to a fractal medium ( Loverdo et al 1. l2009d) . In this case, using 
the lO'Shaughnessv and Procaccial ( 1985() formalism, the large x behavior of the relocation distribution is obtained as 

7r(a;|ro) ^ r^™ j^i+d ^-dj ^ ^j^ere ^1 is the dimension of the projection of the fractal on the axis parallel to the 
cylinder (jLoverdo et al.l . [2009c( ). Hence, the relocation distribution always decays faster than in the case of normal 
diffusion (~ \ / zIt? [z / a)) . As a consequence, the proportion of jumps in the case of random conformation of DNA 
in the fractal type crowding scales like C(^) ^ ^-d^'+df . much smaller than for regular diffusion, which shows 
that fractal crowding favors hops and changes the overall intermittent search. 



d. Beyond the mean : variability of the search time 

Min imizing the mean search time is the optimization procedure the most often used (see for example ICoppev et al.l 
(I2004D 1. However, the entire distribution of the search time can be needed to assess the search kinetics on all time 



scales. Obviously in the case of simple exponential distributions of the search time, the mean is sufficient to describe 
the full dynamics. However, several models discussed below have shown that the search time distribution could not 
always be a single exponential. Such departure from an exponential distribution can have important consequences, 
such as large fluctuations of the search time which could be an extra source of variability in gene expression. 

Effect of trapping sequences 

A first possible source of fluctuations in the search time could come from the existence of trapping sequences along 
the DNA. Since a target sequence is typically 10 base pairs long, very similar sequences are statistically unavoid able and 
can be expected to be local minima in the protein/DNA interaction energy landscape. iBenichou et al.l ( 20091 ) propose 



a model in which the protein can be stuck on such sequences which are similar to the target. The corresponding 
trapping times naturally introduce new time scales in the problem, which potentially could be very long. In the 
framework of this model, it is shown that the search time distribution is best described by two exponentials. In 
particular, the mean search time is controlled by long trapping events even when they are very unlikely, and can be 
orders of magnitude larger than the median search time which is controlled by the trajectories that do not fall into 
the traps. The main outcome of such model is that it reconciles the possibility of having long lived stable complexes 
{ie deep traps), and very fast typical search times, which were casted so far as paradoxical requirements. 

Effect of n searchers 

The influence o f the number o f searc hers has been discussed by s everal authors l Eliazar et al. ] (|2007l iooi); 



iMeroz et aP (|2009t ): ISokolov etall (j^ol). Importantly, as stressed in iBenichou et al.l ([2009t ). the mean time for 



n independent searchers is simply given by the mean time divided by the number of searchers only if the search time 
distribution for a single searcher is a single exponential. In the case of non exponential distributions {eg a sum of 
weighted exponentials) it can be shown that the effect of the number of searchers can be much stronger since the 
weight of long ti me scales decreases ex ponentially with n, selecting for n large enough only the shorter time scale of 
the problem (see IBenichou et aL l (|2009l) 1. Moreover, if the concentration of searchers increases considerably, searchers 



cannot be considered as independent anymore and will act as "roadblocks" along the DNA. These roadblocks will 
decrease the effective sliding length, and may also hide the target, overall slowing down the search process. This 
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results in a trade-off, as stressed in iLi et al. ] (I2009D . On the one hand, the more the proteins, the more the searchers 



for the target and the quicker the search. On the other hand, the more the proteins, the more the crowding, the less 
efficient the search of a single protein. These authors predict that the optimum is obtained for 10"* — 10^ DNA binding 
proteins for E.Coli, which is close to the experimental value of 30 000 proteins. 

Dependence on the starting point: colocalization 

Another origin of fluctuations of th e search time can be due to its dependence on the starting position of the 



protein. IWundcrlich an d MiruYl ( 20081) show that a target which is close to the starting point of the protein can be 



found within a single sliding phase, which yields a very short search time and a rather low variability. In contrast, 
if the target is far away from the starting point of the protein, it is found after numerous 3D excursions. The mea n 
search time is then much long er, and the spread of the distribution of the search time is larger. iKolesov et al.l ( 2007t l: 



IWunderlich and Mirnvl (j2008f) then argue that for increasing the efficiency of the transcr iption factor, it s codin g 
sequence {ie its starting posit ion) should be colocalized with its target sequence (see also iBenichou et al.l ( 2008bf ): 
iBenichou and Voituried (|2009[ ) for a further optimization of this colocalization effect with respect to the diffusion 
coefficient of the protein). Colocalization is indeed observed in real prokaryotes genomes. Such mechanism can 
however be invoked only in prokaryotes, where there is no cell compartments separating protein production from 

DNA. 

More generally, geometric effects on search kinetics have been discussed in ( Benichou et al.l . l2010t ). where it is shown 



that low dimensional effects, such as sliding or diffusion on fractals, can lead to non exponential distributions which 
depend strongly on the starting position of the searcher. Such mechanism could be important for eukaryots. Indeed, 
in eukaryotes, DNA is packed inside the nucleus in what is called the chromatin, and some DNA regions could be 
more or less accessible depending on the chromatin configuration. For example the DNA close to the nucleus pores is 
much more accessible t o incoming protein s than the DNA buried deep inside the nucleus, leading potentially to very 



mucn more accessible t o incoming protein s tnan tne UiNA buried aeep mside tne nucleus, leading potentially to very 
different search times. iKampmannI ( 20051 ) argues qualitatively that proteins binding to DNA could take advantage 



of the heterogeneities of the chromatin, and, depending on the searched sequence, adopt different optimal strategies. 
The geometric effect can be particularly important in the case of genes which need to be activated simultaneously. 
Indeed, their colocalization in the nucleus permits to share the transcription material, since the search time for a 
transcription factor going from one to the other will be much smaller than in the case of a random localization in the 
nucleus. 



4. Conclusion on protein/DNA interactions 

The mechanism of facilitated diffusion of proteins on DNA is intrinsically intermittent in the general meaning 
defined in this review : it is a combination of one dimensional motion in close interaction with DNA, called "sliding", 
which enables target detection, and faster 3D excursions. From the theoretical point of view, this further example of 
intermittent process has been shown to significantly speed up the search. Over the past few years, strong experimental 
evidences have been brought in, showing that this mechanism is indeed at play. Interestingly, in this case in the "fast" 
phase the searcher is not able to detect the target not because of a lowering of its perception abilities, but simply 
because the motion takes place in a geometrical space which does not contain the target. Most of microscopic 
realizations of intermittence fall in this case as we will see in next section on a further example. 



B. Active transport of vesicles in cells 



After this first microscopic example of intermittent search, we turn to another example : active transport of vesicles 
reacting at specific locations in cells. 



1. Active transport in cells 



Various motor proteins such as kinesins or myosins are able to convert the chemical fuel provided by ATP into 
mechanical work by interac ting with the semiflexible oriented filaments (mainly F-actin and microtubules) of the 
cytoskeleton ([Albertsl . [200l . As many molecules or larger cellular organelles like vesicles, lysosomes or mitochondria, 
hereafter referred to as tracer particles, can randomly bind and unbind to motors, the overall transport of a tracer in 
the cell can be described as alternating phases of sta ndard dif fusive tra nsport (sometim es subdiffusive) , and phases of 
active directed tra nsport powered by motor proteins ( Albertsl . [2002; Sal man et aLl , l2005D (see figure [H)) . In particular, 
iHuet et all (|2006l) studied the rate of transitions between ballistic, diffusive and "on the target" states of vesicles, and 
found that the vesicles studied are much more likely to react in the free diffusive phase than when bound to motors. 
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Rigid polymer sucii as 
microtubule or actin 




Does not need ATP 
Distance^ oc Dt (or less if subdifltusive) 
Diffusion coefficient decreases when the size of the object 
increases 



Consumes ATP 
Distance = Vt 
Velocity ~ 1 (j.m.s^^, almos t independent of th e object 
size (see for example iBlock et all (120031 )') 



Figure 14 Transport options for vesicles inside cells. 



Active transport in this case is therefore clearly a further example of intermittent behavior. Active transport in 
cells has been extensively studied both experimentally, for i nstance by single particle tracking methods ( Caspi et al.L 
I2OOOI . I2OO2I : [Howard et aU , 'l989' 'Sheetz and Spudich', 198 3|), and theoretically by evaluat ing the mean displacement 
of a tracer (lAidari. 1995,: Salman ct al. , 2005; .Shlcsinger and Klafted . Il989l ). or stationary concentration profiles 
( Nedelec et al.l , 200 ll ). This transport is important for example for dynamically regulating the distribution of proteins 
such as membrane receptors. 

Most of cell functions are regulated by coordinated chemical reactions which involve low concentrations of reactants 
(such as ribosomes or vesicles carrying targeted proteins), and which are therefore lim ited by t ransport. An analytica l 
model based on the idea of intermittence has been introduced in lLoverdo et all ( 2008[ ) (see also lLoverdo et aL j ((20093); 
iMirnvl (|2008[ )). and enables the determination of the kinetic constant of transport limited reactions in active media, 
and further shows that the kinetic constant can be optimized. We give below the main results of the model, and 
further details can be found in appendix. 



2. Model 

The model relies on the idea of intermittent search strategies and has important similarities with the model of 
section HH which are discussed in section IIVI We consider a tracer particle evolving in a d-dimensional space (in 
practice d = 1,2,3) which performs thermal diffusion phases of diffusion coefficient D (denoted phases 1), randomly 
interrupted by ballistic excursions bound to motors (referred to as phases 2) of constant velocity V and direction 
pointing in the solid angle wv (see figure [TS)) . The distribution of the filaments orientation is denoted by p{uj\-), and 
will be taken as either disorde red o r pola rized (see figures [121 [13 fT9)) , which schematically reproduces the different 
states of the cytoskeleton (Albertsl . |2002[ ). The random duration of each phase i is assumed to be exponentially 
distributed with mean r^. The tracer T can react with reactants R (supposed immobile) during free diffusion phases 1 
only, as T is assumed to be inactive when bound to motors. Reaction occurs with a finite probability per unit of time 
k when the tracer-reactant distance is smaller than a given reaction radius a. In what follows the kinetic constant K 
of the reaction T + R ^ R is explicitly determined. 



3. Methods 

We now present the basic equations in the case of a reactant centered in a spherical domain of radius b with reflecting 
boundary. This geometry both mimics the relevant situation of a single target and provides a mean field approximation 
of the general case of randomly located reactants with concentration c — a'^/b'^, where b is t he typical distance b e tween 
reactants. We start from a mean field approximation of the first order reaction constant ( Berg and Blomberd . Il976[ ) 



and write K = 1/ ( t), where the key quanti t y of our appro ach is the reaction time (t) which is defined as the mean 
first passage time (jCondamin et al. I l2007bl : iRed^ l200lh of the tracer at a reactant position uniformly averaged 
over its initial position, ti is defined as the mean reaction time if the tracer starts in phase 1 at position r, and 
t2 is defined as the mean reaction time if the tracer starts in phase 2 at position r with velocity v. F or the active 
intermittent dynamics defined above, ti and <2 satisfy the following backward equations ( Redneii [200 1[) (see section 
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IIVI for derivation) : 



DA„ti + — (t2- ti)p{ujv)dLu^r - kla{r)h = -1 



Tl „ 

V-Vr^a {t2-ti) = -l 

T2 



(36) 



where Ar and Vr are the Laplacian and the gradient on the initial position, is the indicator function of the sphere 
of radius a. As these equations p6|) are of integro-differential type, standard methods of resolution are not available 
for a general distribution p. 

However, in the case of a disordered distribution of filaments {pioJv) = ^/^d-, where Via is the solid angle of the 
d-dimensional sphere), these equations can be solved exactly in dimension 1. In dimension 2 and 3, an approximate 
scheme has to be introduced; the details of the calculation are given in appendix. We present here simplified expressions 
of the resulting kinetic constant by taking alternatively the limit — >■ oo, which corresponds to the ideal case of perfect 
reaction, and the limit Z? -> which allows us to isolate the k dependence. 



4. Active transport in the cytoplasm 



Unbinding 




Figure 15 Vesicle transport in the bulk (3 dimensions). 



We first discuss the d = 3 disordered case (see figure [T5|). which provides a general description of the actin cytoskele- 
ton of a cell in non polarized conditions, or of a generic in vitro active solution. An analytical form of the mean first 
passage time (t) =1/ Kj,d is given in appendix, and plotted in figure [161 Strikingly, K^d can be maximized as soon as 
the reaction radius exceeds a threshold Cc — QD /V for the following value of the mean interaction time with motors: 



opt 

^2, 3d ~ 



•\/3a 



1.078f 



(37) 



where xq is the solution of 2tanh(a;) — 2x + a;tanh(a;)^ = 0. The ti dependence is very weak, but one can roughly 



estimate the optimal value by t°3*^ — &D/V^. This gives in turn the maximal reaction rate 



^3d — 



cV \/i (xq — tanh(a;o)) 



(38) 



so that the gain with respect to the reaction rate iff^ in a passive medium is G^d = K^/K^j^ ~ CaV/D with C ~ 0.26. 

Several comments are in order, (i) First, 3*^ neither depends on D, nor on the reactant concentration. A 
similar analysis for k finite (in the D ^ limit) shows that this optimal value does not depend on k either (see 
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b/a=5 b/a=40 



Figure 16 Optimization of the reaction rate for intermittent active transport. Gain of reactivity due to active transport in 3 
dimensions as a function of T2 for different values of tiie ratio b/a (logarithmic scale). The analytical form (the mean detection 
time with diffusion alone divided by the mean detection time with intermittence) (plain lines) is plotted against numerical 
simulations (symbols) for the following values of the parameters (arbitrary units): a — 1 (green, □), a = 5 (blue, *), a = 7 
(purple, o), a — 10 (red, +), a — 14 (brown, x), a = 20 (orange, o), with ti = 6, V = 1, D = 1. G-id presents a maximum only 
for a > Qc — 4. 



section irv|) . which proves that the optimal mean interaction time with motors is widely independent of the parameters 
characterizing the diffusion phase 1. (ii) Second, the value Oc should be discussed. In standard cellular conditions 
D ranges from ~ 10^^ [a.m^.s~^ for vesicles to ~ 10 |j,m^.s~^ for small proteins, whereas th e typical veloci ty of a 
motor protein is ~ 1 la^m.s^^, value which is widely independent of the size of the cargo (|Albertd [200^ . This 
gives a critical reaction radius Oc ranging from ~ 10 nm for vesicles, which is smaller than any cellular organelle, 
to ~ 10 |a.m for single molecules, which is comparable to the whole cell dimension. Hence, this shows that in such 
3-dimensional disordered case, active transport can optimize reactivity for sufficiently large tracers like vesicles, as 
motor mediated motion permits a fast relocation to unexplored regions, whereas it is inefficient for standard molecular 
reaction kinetics, mainly because at the cell scale molecular free diffusion is faster than motor mediated motion. This 
could help justifying that many molecular species in cells are transported in vesicles. Interestingly, in standard 
cellular conditions t^"^^^ is of order 0.1 s for a typical reaction radius of order 0.1 jam. This value is compatible with 
experimental observations ( Albertsl . |2002| ). and suggests that cellular transport is close to optimum, (iii) Last, the 
typical gain for a vesicle of reaction radius a > 0.1 [a.m in standard cellular condit ions is G^d ^, 2 . 5 (see figure [T51) and 
can re ach G^d — 10 for the fastest types of molecular motors (1/ ~ 4 |j.m.s~^, see I Albertsl (|2002D : ISheetz and Soudichl 
( 19831 )). independently of the reactant concentration c. As we shall see below the gain will be significantly higher in 
lower dimensional structures such as axons. 



5. Active transport at membranes 

We now come to the d = 2 disordered case (see figure [T7)) . Striking examples in cells are given by the cytoplasmic 
membrane, wh ich is closely coupled to the network of cortical actin filaments, or the lamellipodium of adhering cells 
(I Albertsl . [200^ . n many cases the orientation of filaments can be assumed to be random. It can be shown that as 
for d = 3 (see section WV} . the reaction rate K2d can be optimized in the regime D/V <C a <C 6. Remarkably, the 
optimal interaction time T2 2*^; takes one and the same value in the two limits fc — >■ cx) and -D — >■ : 



(39) 
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Figure 17 Planar structures such as membranes and lamellipodia (d — 2). 



which indicates that again 2*,; does not depend on the parameters of the thermal diffusion phase, neither through 



D nor k. In the limit fc — > oo one has r-^ 



opt 
2d 



in"i/c)-i ' ^'^^ maximal reaction rate can then be obtained 



k: 



cV 



2d 



av/21n(l/c) 



(40) 



Comparing this expression to the case of passive transport yields a gain G2d = ^m/^m — o,v ^y\n{l / c) / {4D^/2) . 
As in the c? = 3 case, this proves that active transport enhances reactivity for large enough tracers (with a critical 
reaction radius ac ~ D/V of the same order as in the d ~ 3 case) such as vesicles. However, here the gain G2d depends 
on the reactant concentration c, and can be more significant : with the same values of D, V and a as given above 
in standard cellular conditions, and for low concentrated reactants (like specific membrane receptors) with a typical 
distance between reactants b> 10 \xm, the typical gain is G2d 8, and reaches 10 for single reactants (like examples 
of signaling molecules) (see figure [TS]) . 



6. Active transport in tubular structures 



The case of nematic order of the cytoskeletal filaments, which depicts for instance the situation of a polarized cell 
(lAlbertsl . [200^ . can be shown to be equivalent in a first approximation to the 1-dimensional case, which is exactly 
solvable (see figure [T^ (for calculations, see section lA.l.b[) . The d = 1 case is also important on its own in cel l biology 
as m any 1-dimensional active structures such as axons, dendrites, or stress fibers are present in living cells ( Alberta . 
I2OO2I ). As an illustration, we take the example of an axon, filled with parallel microtubules pointing their plus end 
in a direction e. We consider a tracer particle interacting with both kinesins ("+" end directed motors, of average 
velocity Ve ) and dyneins ("-" end directed motors, of average velocity —Ve) with the same characteristic interaction 
time T2 (see figure lb). For this type of tracer, the mean first passage time satisfies equations ([55]) with an effective 
nematic distribution of filaments p{uJv) = ^{S(V — e) + (^(V + e)). The reaction rate Kij, is maximized in the regime 
D/V <^ a b for the following values of the characteristic times (see figure [^0)1 : 



opt 
'l,ld 



1 D 

48 V^' 



opt 
^2.1d 



V3^ci/2' 



for fc — > 00. The maximal reaction rate KVl is then given by 



K 



Id 



2a 



(41) 



(42) 



and the gain is Gid = K^JK^^ ~ aV/{2^/SDc^^'^), which proves that active transport can optimize reactivity as 
in higher dimensions. Very interestingly the c dependence of the gain is much more important than for d = 2,3, 
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Figure 18 Optimization of the reaction rate for intermittent active transport. Gain of reactivity due to active transport G2d in 
two dimensions as a function of ri or T2 (logarithmic scale). The analytical form (the mean detection time with diffusion alone 
divided by the mean detection time with intermittence (|A50|) ) (plain lines) is plotted against numerical simulations (symbols) 
for the following values of the parameters (arbitrary units): a = 20, b — 2000 (brown, a — 10, b = 1000 (red, □), a = 10, 
b = 100 (green, +), a — 2.5, b — 250 (blue, o) with V = 1, D = 1. These curves represent standard cellular conditions (as 
discussed in the text). 




Figure 19 Tubular structures in cells such as axons and dendrites {d — 1). 



which shows that the efficiency of active transport is strongly enhanced in l-dimensional or nematic structures at low 
concentration. Indeed, with the same values of D, V and a as given above in standard cellular conditions, and for a 
typical distance between reactants b > 100 |i.m (like low concentrated axonal receptors), one obtains a typical gain 

Gid ^ 100 (see figure [20l) . In the limit of finite reactivity (fc finite and D ^ 0) one has r^*^'^ ~ Vvk and 
the same optimal value (j4l]) of T2^id- As in higher dimensions '''2'^*^ depends neither on the thermal diffusion coefficient 
D of phases 1, nor on the association constant k, which shows that the optimal interaction time with motors t!^^^ 
presents remarkable universal features. Furthermore, this approach permits an estimate of Tj''' compatible with 
observations in standard cellular conditions, which suggests that cellular transport could be close to optimum. 
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Figure 20 Optimization of tlie reaction rate for intermittent active transport. Gain of reactivity due to active transport Gid in 
one dimension as a function of ri or T2 (logaritlimic scale). Tiie analytical form (the mean detection time with diffusion alone 
divided by the mean detection time with intermittence (|A23|) ) is plotted against the exact solution (symbols), for the following 
values of the parameters (arbitrary units): D = 1,V = 1 for all curves and a = 10, 6 = 10"* (red, +), a = 10, b = 10^ (blue, o), 
a = 2.5, b = 10'^ (green, □). Standard cellular conditions (as discussed in the text) correspond to blue and red curves. 



7. Conclusion on intermittent active transport 

Starting from the observation of vesicles alternating free diflFusion and phases bound to motors performing ballistic 
motion, and from the observation that (at least in some cases), vesicles can only react in the free phase, a model for 
intermittent active transport has been proposed. The reaction rate, which can be approximated by the inverse of the 
mean first passage time, can be explicitly calculated in this model for various cellular geometries (bulk cytoplasm, 
membranes, tubular structures). This shows that intermittent transport can indeed increase reaction rates, in partic- 
ular for large objects such as vesicles, and in particular in low dimensions. The model for the reactive phase is either 
diffusive or static (with a reaction rate) , and both lead to the same optimal duration of the ballistic phase. The latter 
point is investigated in more details in next section. 



IV. INTERMITTENT SEARCH : A ROBUST STRATEGY 

As shown in previous sections, intermittent search strategies are observed at the macroscopic scale (foraging animals) 
as well as at the microscopic scale (localization of a DNA sequence by a protein, vesicle transport in cells) . The models 
we h ave used to interpret these findings, in particular in sections III. CI and IIII.Bl present similar general features. 

In lLoverdo et aL l (l2009bD . a generic model of intermittent search based on these general features has been introduced, 
and studied systematically in 1, 2, and 3 dimensions, and for three different modelings of the detection phase. This 
rather technical section (completed by appendix) gathers the main tools usually involved in the calculation of first- 
passage properties of intermittent random-walks, and utilized throughout this review. Finally, general conclusions on 
intermittent random walks can be drawn from this systematic study, and are summarized in table HI 



A. Introduction 

The generic model presented in this section follows the original definition of intermittence given in introduction 
and relies on a succession of slow phases with detection, and ballistic phases without detection, without direction 
correlation between ballistic phases. This model is minimal in the sense that the searcher has low memory skills. 
Indeed, without correlations between ballistic phases, there is no spatial memory. We also assume a Markovian 
searcher, i.e. with no temporal memory. As previously we address the following main questions: is it beneficial for 
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the search to perform such fast but non reactive phases? Is it possible, by properly tuning the kinetic parameters of 
trajectories (such as the durations of each of the two phases) to minimize the search time? We develop in what follows 
a systematic analytical study of intermittent random walks in one, two and three dimensions and fully characterize 
the optimal regimes. Overall, this systematic approach allows us to identify robust features of intermittent search 
strategies. In particular, the slow phase that enables detection is often hard to characterize experimentally. Here we 
propose and study three distinct modelings for this phase, which allows us to assess to which extent our results are 
robust and model independent. Our analysis covers in details intermittent search problems in one, two and three 
dimensions and is aimed at giving a quantitative basis - as complete as possible - to model real search problems 
involving intermittent searchers. 

We first define the model and introduce the methods. Then we summarize the results for the search problem in 
dimension one, two and three, for different types of motion in the slow phase. Eventually we synthesize the results in 
the table H] where all cases, their differences and similarities are gathered. This table finally leads us to draw general 
conclusions. 



B. Model and notations 

1. Model 



We consider an intermittent searcher that switches between two phases. The switching rate Ai (resp. A2) from 
phase 1 to phase 2 (resp. from phase 2 to phase 1) is time-independent, which assumes that the searcher has no 
temporal memory and implies an exponential distribution of durations of each phase i of mean = l/A^. 



phase 1 



phase 1 



"■-^ phase 2 



phase 2 





Static mode 



Diffusive mode 




Ballistic mode 



Figure 21 The three different descriptions of phase 1 (the phase with detection), here represented in two dimensions. 



Phase 1 denotes the phase of slow motion, during which the target can be detected if it lies within a distance 
from the searcher which is smaller than a given detection radius a, which is the maximum distance within which the 
searcher can get information about target location. We propose 3 different modelings of this phase, in order to cover 
various real life situations (see figure [?T|) . 

• In the "static mode", the searcher is immobile, and detects the target with probability per unit time k if it lies 
at a distance less than a. 

• In the second modeling, called the "diffusive mode", the searcher performs a continuous diffusive motion, with 
diffusion coefficient D, and finds immediately the target if it lies at a distance less than a. 

• In the last modeling, called the "ballistic mode", the searcher moves ballistically in a random direction with 
constant speed vi and reacts immediately with the target if it lies at a distance less than a. W e note that this 
mode is equivalent to the model of Levy walks searches proposed in IViswanathan et al.l ( 1999h . except for the 
law of the time between reorientations (see section III.A|) . It was shown that for destructive search, i. e. targets 
that cannot be revisited, the optimal strategy is obtained for a straight ballistic motion, without reorientations 
(see section Hi. Ap . In what follows it is shown that if another motion, "blind" {i.e. without detection) but with 
higher velocity is available, there are regimes outperforming the straight line strategy. 

Some comments on these different modelings of the slow phase I are to be made . First, these three modes schematically 
cover experimental observations of the behavior of animals searching for food ( Belli . [iQQlllO'Brien et ah I I1990D . where 
the slow phases of detection are often described as static, random or with slow velocity. Several real situations 
are also likely to involve a combination of two modes. For instance the motion of a reactive particle in a cell not 
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bound to motors can be described by a combination of the diffusive and static modes. For the sake of simphcity, 
these modes are treated independently, and this approach can therefore be considered as a Hmit of more reahstic 
models. Finally, combining these three schematic modes covers a wide range of possible motions, from subdiffusive 
(even static), diffusive, to superdiffusive (even ballistic). Beyond the modeling of real- life systems, studying different 
detection modes enables us to assess the robustness of the results. 

The phase 2 denotes the fast phase during which the target cannot be found. In this phase, the searcher performs a 
ballistic motion at a constant speed V in a random direction, redrawn for each new phase 2, independently of previous 
phases. In real examp les correlations between successive ballistic phases could exist, as observed for foraging animals 
(lO'Brien et allllOQd) . ' f correlations are very high, it is close to a one dimensional problem with all the phases 2 in 
the same direction, a different problem already treated in section III.Bl We consider here the limit of low correlation, 
that is of a searcher with no memory skills. 

We assume that the searcher evolves in a d-dimensional spherical domain of radius 6, with reflective boundaries 
and with one centered immobile target of radius a. As the searcher does not initially know the target's location, 
we start the walk from a random point of the d-dimensional sphere, and average the mean target detection time 
over the initial position. This geometry models the case of a single target in a finite domain, and also provides a 
good approximation of an infinite space with regularly spaced targets. Such regular array of targets corresponds 
to a mean-field approximation of random distributions of targets, which can be more realistic in some experimental 
situations. 



2. Methods 

We explain here the general methods, and introduce the notations. 

We define Si{r,t) the probability that the searcher has not yet found the target at t, starting from r in state i, 
where state i = 1 is the slow motion phase with detection and state i = 2 is the fast motion phase without target 
detection. Note that in dimension one, the space coordinate will be denoted by x, and in the case of a ballistic mode 
for phase I, the upper index in tf stands for ballistic motion with direction ±x. The su rvival probabilit y Si(r,t) is 
solution of the following standard backward differential Chapman- Kolmogorov equations (jGardinerl . Il996[ ) : 

Lls,ir,t) + - {s,{r,t) - s,{r,t)) - ks,{r,t)Ia{r)S{i - 1) = (43) 
n at 

with /a(r) = 1 when r < a, and else, and l\ is the adjoint operator of the transport operator. For example l\ = DA 
for diffusion and i| = v • V for a ballistic motion of velocity v . The mean first passage time to the target i,; (r) for a 
searcher starting in the phase i from point r is then given by : 

t,{v) = - t^^^^dt = s,{v,t)dt (44) 

Consequently, for each phase i, ti{r) is solution of : 

L|t,(r) + - (tj(r) - t,(r)) - kt,{r)Ia{r)S{i - 1) = -1. (45) 

We assume that the searcher starts in phase 1, and to take into account the fact that it does not initially know the 
target's location, we average the mean detection time over the starting point, leading to the following definition of 
the mean search time : 

^'" = "i77?rT / ^iWrfr, (46) 

with Qd the d-dimensional sphere of radius b and V{fld) its volume. Unless specified, we will consider the low target 
density limit a <C &. 

Our general aim is to minimize tm as a function of the mean durations ti , T2 of each phase, and in particular 
to determine under which conditions an intermittent strategy (with finite T2) is faster than a usual single state 
search in the phase 1 only, which is given by the limit ri — > 00. In the static mode, intermittence is necessary for 
the searcher to move, and is therefore always favorable. In the diffusive mode, we will compare the mean search 
time with intermittence tm to the mean search time for a single state diffusive searcher idiff, and define the gain as 
gain = fdiff/^m- Similarly in the ballistic mode, we will compare to the mean search time for a single state ballistic 
searcher t^ai and define the gain as gain = tbai/tm- 
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The upper index "opt" is used to denote the value of a parameter or variable at the minimum of tm . 

Calculations are exact in dimension one, whereas approximation schemes (which can be checked numerically) are 
needed in dimen sion two and three. T he main calculation steps are given in appendix and further technical details 
can be found in (|Loverdo et all . l2009bl ) . We now summarize the main results for each dimension. 



C. Dimension one 



Besides the fact that it involves more tractable calculations, the one-dimensional case is also interesting to model 
real search problems. As discussed before, at the micros copic scale, tu bular structures of cells such as axons or 
dendrites in neurons can be considered as one-dimensional ( Albert si . l2002l ). The active transport of reactive particles, 
which alternate between diffusion phases and ballisti c phases when bound to molecular motors, can be schematically 
captured by this generic model with diff usive mode ( Loverdo et al.l , [20081 ). At the macroscopic scale, one could cite 
animals like ants ( Dussutour et al.l , l2005l ) which tend to follow tracks or one-dimensional boimdaries. More generally, 
borderlines between different habitats, such as a shoreline, can be considered as one-dimensional. 

It is showed in appendix that intermittent search strategies in dimension one share similar features for the static, 
diffusive and ballistic detection modes. In particular, all modes show regimes where intermittence is favorable and 
lead to a minimization of the search time. Strikingly, the optimal duration of the non-reactive relocation phase 2 is 

quite independent of the modeling of the reactive phase : r2 ^* = for the static mode, for the ballistic mode 



(in the regime vi < vf ~ "j y x)' ^-'^'^ ^^"^ diffusive mode (in the regime b > y and a ^ V'\J a)' This shows the 
robustness of the optimal value ^* . 



D. Dimension two 



As discussed previously, the two-dimensional problem is particularly well suited to model animal behavi ors; it is 
also r elevant to the microscopic scale, since it mimics for example the case of cellular traffic on membranes (| Albertsl . 
l2002f ). While in dimension one the mean search time can be calculated analytically, in dimension two (and later in 
dimension three) approximation schemes are necessary, and can be checked by numerical simulations. 

Remarkably, for the three different modes of detection (static, diffusive and ballistic), there is a regime where inter- 
mittence minimizes the search time for one and the same ^* , given by — y y^ln (^) — i. As in dimension one, 

this indicates that optimal intermittent strategies are robust and widely independent of the details of the description 
of the detection mechanism. 



E. Dimension three 



The three dimensional case is also relevant to biology. At the microscopic scale, it corresponds for example to 
intracellul ar traffic in the bulk cy toplasm of cells, or at larger scales to ani mals living in dimension thr ee, such as 
plankton ( Bartumeus et all , l2003f ) . or C.elegans in its natural habitat (soil) (jKiontke and Sudhausl . [2b05[ ) . As it was 
the case in dimension two, different assumptions have to be made to obtain analytical expressions of the search time. 
Such assumptions can be checked by numerical simulations using the same algorithms as in dimension two. 

For the three possible modelings of the detection mode (static, diffusive and ballistic) in dimension three, there is 
a regime where the optimal strategy is intermittent. Remarkably, and as was the case in dimensions one and two, 
the optimal time to spend in the fast non-reactive phase 2 is independent of the modeling of the detection mode and 
reads r2^* ~ 1-1 "p^- Additionally, while the mean first passage time to the target scales as 6^, the optimal values of 
the durations of the two phases do not depend on the target density a/b. 



F. Discussion and conclusion 

To summarize, the methods of calculation developed in this section allow one to show that the mean search time 
of intermittent random walks can be minimized under broad conditions. The table |l] summarizes the main results 
of this minimization. This study shows that the optimal durations of the two phases and the gain of intermittent 
search (as compared to a single state search) do depend on the target density in dimension one. In particular, the 
gain can be very high at low target concentration. Interestingly, this dependence is smaller in dimension two, and 
vanishes in dimension three. The fact that intermittent search is more advantageous in low dimensions (1 and 2) can 
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be understood as follows. At large scale, the intermittent searcher effectively performs a random walk, and therefore 
scans a space of dimension two. In an environment of dimension one (and marginally of dimension two), the searcher 
therefore oversamples the space, and it is favorable to perform large jumps to go to previously unexplored areas. On 
the contrary, in dimension three, the random walk is transient, and the searcher on average always scans previously 
unexplored areas, which makes large jumps less beneficial. 

Additionally, these results show that, for various modeling choices of the slow reactive phase, there is one and 
the same optimal duration of the fast non reactive phase, which depends only on the space dimension. This further 
supports the robustness of optimal intermittent search strategies. Such robustness and efficiency could explain why 
intermittent trajectories are observed so often, and in various forms. 
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Table I Recapitulation of main results of the generic intermittent search model : strategies minimizing the mean first passage 
time on the target. In each cell, validity of the regime, optimal ri, optimal T2, minimal tm (tm with rt = t°^*). Yellow 
background highlight the value of Tj*'' independent from the description of the slow detection phase 1. Results are given in the 
limit 6 3> a. Complements to the intermediate regime for the diffusive mode in 2 dimensions : c = 4 (7 — ln(2)) with 7 the Euler 
constant; w solution of In (41n(w) - 5 + c) = -8(lnw)^ + (6 + 8 ln(fe/a)) ln(«;) - 101n(fe/a) + 11 -c(c/2 + 21n(a/6) - 3/2); 

[1 + 2 \n{w) In (^)] . Adapted from lLoverdo et all (|2009bl ') 



in this regime we have t° 



■In I 



a / 4 \n{w) — 5 



1 + 



mg(41n(m)-7) 
hV ^ 5 In(ui) —5 
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V. EXTENSIONS AND PERSPECTIVES 

Far from closing the problem, the generic model presented previously opens interesting perspectives. In this section, 
we highlight a few promising directions : (i) Influence of the targets distribution ; (ii) Taking into account a more 
involved searcher, enjoying now some orientational and temporal memory. Indeed, in the generic model of section 
IIVI the searcher has minimal memory skills. On the one hand, the phase duration distribution is exponential, which 
means that there is no temporal memory : the effect of other duration distributions is studied in section IV. CI On 
the other hand, the direction of each new ballistic "blind" phase is taken at random, independently of the previous 
phases, meaning that there is no orientational memory : we study the effect of correlations in section IV. Bl (iii) 
Moving targets, which can be more realistic both at the microscopic and macroscopic scale. Next, we briefly review 
in section IV. El similar models of intermittent search which have been proposed recently in other contexts and finally 
discuss how further models could also be applied to design efficient searches instead of interpreting biological systems 
(see section IV.FI) . 



A. Influence of tfie target distribution on the search time 

We first study the influence of target distribution on the previous results. For the sake of simplicity we study the 
one dimensional model of section III.BI 



1. How real targets are distributed? 



In the context of foraging animals, target distributions are often described as regular, random or patched ( Belli . 



I1991D (see fi gure . In the models presented in previous sections, the chosen geometry can be interpreted as one 
target in a finite domain, or as an infinite array of regularly spaced targets. The regular distribution is representative 
of the real-life case of targets that repel each other, thus being as far from each other as possible. This distribution 
is also a mean-field approximation of other distributions. As the regular distribution has already been studied, let us 
discuss the other representative distributions. 



regular distribution 



• — • ■ — • ■ — — 



Poissonian distribution 

mm mm m-mm 

patched distribution 

Figure 22 Examples of target distributions. 



If targets are in patches, for example when they attract each other, when a target is found it is likely that other 
targets are present in the immediate sur roundings. Thus a simple strategy is to switch behavior when a target is 
encountered, as proposed for example bv iBenhamorJ ( 19921 ). The search is then in two steps : finding a patch, and 
exploiting it. For the first step, previous results are still valid, except for the density of targets which has to be 
replaced by the density of patches. 

In the following we focus on the last case of Poissonian targets, which corresponds to situations of non-interacting 
targets. 



2. Analytical results in the case of a Poissonian distribution of targets 

In the case of a 1 dimensional Poissonian distribution of targets, the distance between two consecutive targets is 
exponentially distributed. Except for this change, the other parameters remain as defined in the model of section 
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The mean search time is in general hard to calculate for a Poissonian target dis tribution, whi c h can be see n as 
frozen disorder. However, estimates (for L ^ D/V) can be given in 3 regimes fsee iMoreau et al ] ([2Q03, [20093) for 
details) : 

• In the large ballistic displacements limit (when Vt2 ^ ^/Dti) , two successive diffusive phases can be considered 
as non-overlapping. It can be shown that in this regime : 

(47) 



• In the small ballistic displacements limit (when Vt2 <C V Dti ) , successive diffusive phases often overlap. It leads 
to : 

W-i^. (48) 

• The most interesting situation is the intermediary regime. Indeed, in the first case (large ballistic displacements), 
relocations are too long and overshoot the target; and in the second case (small ballistic displacements), there 
are often repetitive scans of the same areas. In the intermediary regime, the mean first passage time to the 
target can be approximated by : 

n+r^ (1+^)^(1 + 6^ 
- ^ Vt2 (l + 40 + 2e02)' (49) 

with = VT2/y/lM and e = y/Drl/L. 

This last regime enables a discussion of the efficiency of the intermittent search. The efficiency can be quantified by 
comparing (i)"^*, the mean search time with intermittence at the minimum with rjiff = L^/(2D), the mean search 
time with diffusion alone. It can be shown that intermittence decreases the search time in the limit of low target 
density, and that the mean search time is minimized for ti as small as possible. The optimization with respect to 
T2 leads to two regimes, depending on the minimal value of ri as compared to the previously introduced timescale 
T = D jV"^ ^ characteristic of the searcher. : 

• When Tx~> T : 



At the optimum, the mean search time is 



and the gain is 



(50) 



where T-^ai ~ L/V is the typical time needed to travel in the ballistic mode the distance between two consecutive 
targets. As we shall see in the following, in this regime the approximations are very accurate. 



• When Ti ^ T : 



At the optimum, the mean search time is 



opt 
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and the gain is : 



As we shall see in the following, the approximations are qualitatively good in this regime, but not as precise as 
in the other regime. Indeed, the gain obtained here would mean that the mean first passage time to the target 
is smaller than Ti,ai, which is the minimal mean time to travel to the target (except if Ti,ai > T^iff). In fact, as 
can be seen in figure simulations show that {t)°P* -> T^ai. It means that very fast intermittence enables the 
searcher to retain the best of the two phases : reactivity of phase 1 and motion of phase 2. 




T 1-™! '-'-'"T 1-™! T ' I "I 1 I "I ' ' "I 1 I'M ' ' I i| 1 I'M ' ' ' i| 

10"^ 10"3 10"2 10"! lO" 10^ 10^ 10^ 10* 10"^ 10"^ 10"! lO" 10' 10^ 10^ 10* 

ln{{T) / Tdiff) as a function of /n(ri). T2 — 10 (red, □), ln[{T) / r^is) as a function of ln{T2). t\ — 100 (red, □), 
T2 — 1 (green, +), T2 = 0.1 (blue, o). ri = 1 (green, +), n = 0.01 (blue, o). 



Figure 23 Validity of the approximations. Mean first passage time to the target, renormalized by the mean first passage 
time without intermittence. Small ballistic displacements approximation 1)48^ (dashed line). Large ballistic displacements 
approximation (|47p (dotted line). Intermediary approximation (|49p (line). Numerical simulations (symbols). D = 1, V = 1, 
L = 10^ 



Figure [23] represents the mean search time (t) as a function of ri and T2 for typical values of the other parameters. 
It enables to compare the numerical results with the approximations (|47p and (|48p . and with the intermediary 
approximation (|49p . It shows that the approximations of large and small ballistic displacements are valid in the 
expected conditions, and the intermediary approximation (j49p correctly reproduces the existence and the position 
of the minimum of (t) . Figure [24] supports the scaling laws relating ri and the corresponding optimal waiting time 
T2 at the optimum. The exponent 3/4 of the theoretical scaling law (|50p for r <C ri is very well confirmed by the 
simulations. This is not the case for the law ([5^ for t ti, which indicates that the approximations should be 
handled with care for short waiting times ti, T2, although their results are qualitatively correct. Figure [55] shows 
the gain as a function of ri in different possible conditions. It supports the conclusions of the theoretical study, and 
indeed confirms that the gain due to intermittence can be very important if T(jiff 3> r^a;. 



3. Conclusion 



In the case of a Poissonian distribution of targets, intermittence remains valid as a strategy minimizing the search 
time. The optimal strategy still consists in taking ti as small as possible. However, T2^* is different from the 

case of regularly spaced targets. The optimal mean duration of ballistic flights scales as \J^{tit)^^^ in the limit 
Ti ^ T = D/V^. In this regime, at the optimum, (t) ~ hvy t7' "^i^h a gain compared to diffusion alone oc L/ ^/Dt\. 
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Figure 24 ln{T2^^) as a function of ?n(ri). Small ri analytical prediction (|53p (dashed black line). Large n analytical prediction 
pi)) (solid black line). Numerical values (symbols), for L = 10 (green □), L = 10^ (red +), L = lO''^ (blue o). D = 1, V ^ 1. 




Figure 25 ln{G) as a function of /n(Ti) (r2 taken optimal). Small ri analytical prediction (|55p (dotted line). Large n analytical 
prediction H52p (line). Numerical simulations (points). L = 10^ (green, □), L = 10^ (red, +), L = lO^^ (blue, o). D = 1, 1/ = 1. 



B. Taking into account partial correlations in ballistic phases 

1. Motivation 

The models developed in sections III.BI and IIVI are very similar. The searcher alternates between a slow reactive 
phase, and a fast ballistic blind phase. The main difference is that in section III.BI ballistic phases are always in 
the same direction, whereas on the contrary, in section ITVj the direction of each new ballistic phase, is random and 
independent of the previous ballistic phase. 
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In the case of animal trajectories, the successive directions of ballistic phases are usually correlated ( O'Brien et al.l , 
119901 ). We have considered so far two extremes cases : no correlation or infinite range correlations. In both cases, 
there are regimes where intermittence is favorable. However, in the case of infinite range correlations, the shorter the 
duration of each phase, the smaller the search time. In contrast, in the case without any correlation, the minimal 
search time is obtained for finite values of ri and T2, which even diverge with the system size. In the intermediate case 
of finite range correlations, determining the nature of the minimum is an interesting theoretical question. Besides, 
as real biological systems often present correlations, it is an important issue to take into account correlations in 
the generic model of intermittence. We present in what follows the simplest case of the static mode of detection in 
dimension one . 



2. Model 



The searcher is either in the reactive phase 1 (where it is immobile, and finds the target with probability per unit 
time k if the target is at a distance smaller than a), or in the ballistic phase 2, of velocity V. For each new ballistic 
phase, the direction of V is the same as in the previous ballistic phase with probability p, and in the opposite direction 
with probability 1 — p. The distribution of the duration of the phases is exponential, of mean , and the distance 
between two targets is 2b. 

In the case of no correlations {p ~ 1/2) the mean search time has been calculated in the section lA.l.al (see equation 

(jA13p ). where it was shown that the optimum is obtained for r^"^* — ^^Yk (ila)^^^ ^^^"^ '''2^* ~ V\fi.' 

In the general case, the methods of section ITVl can be adapted to calculate analytically the mean search time starting 
from a random position in state 1, which can be written as : 

I b~a ( 2{l~p){b-a)^ 1 u(b-a) ( ^/¥F[ua \\\ 

with u = ^^2(1 — p) + kri. 



3. Minimization of the mean search time 



a. Case of infinite range correlation p = 1 




Figure 26 Minimization of the mean search time for the static mode with infinite correlation (p = 1). a°''* as a function of 
w. Theoretical expression for w small (brown dots), for w intermediate (black solid line), for w large (colored dashed lines). 
Optimization of the exact expression (with r2 — !> and ri = 0x2) (symbols), b — 100 (green, □), & = 10^ (blue, +), b = 10" 
(violet, o), b = 10^ (red, *). a = 1, V ^ 1. 



It can be shown that the mean search time is minimized for ti and T2 tending to 0, with ti — ari- We define 
w = ak/v^ and depending on this parameter : 
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• I <w <w* : a°P* 



2w ' 



• w>w* : a°P* ~ 



w* is defined as the solution of ^ — , — 

minimization of the exact expression of the mean search time (see figure 



These expressions are in good agreement with the numerical 



b. Case of intermediate correlations 

The mean search time obtained in (|56p is difficult to optimize. An important question raised is to determine whether 
the mean search time is minimized for finite ti and T2 (as in the case p ~ 0.5), or for t\ and T2 tending to (as in the 
case p = 1). An answer can be obtained by noticing that a lower bound of the mean search time is given by : 

tm > (n + r,) ^— + ^ 3^.^. ) - (57) 

Supposing that the minimum is realized for at least one of the — > 0, three cases arise. 

• Ti — > with Ti <C T2. In this case i,„ > T2/ (kri) — >■ 00. 

• T2 ^ with T2 ^ Ti . In this case t„i > ti |^^\"^2 00 . 

• Ti ~ T2 and both 0. In this case i,„ > T2 -^^^^p- ~ l/r2 — > 00. 

Finally, this shows that the minimum is realized for finite values of ti and T2 as soon as p < 1. Actually, it can be 
shown that, except for p close to 1, the minimum of the search time is obtained for : 



opt _ I a fb{l-p)^ 



and (58) 



rT = f /P^- (59) 

Interestingly, note that the relation = 2k (r^*^*) obtained initially in the case of the absence of correlations (see 
Table HI still holds in this case. 

These expressions are in agreement with the numerical minimization of the exact expression of the mean search 
time (see figure [57]), except when 1 — p is very small. 



4. Conclusion 



In the simple case of the static mode in one dimension, the influence of correlations on the mean search time and 
its minimization can be studied. An exact expression of the mean search time shows that it is minimized for finite 
values of ri and T2 as soon as p < 1. When {1 — p) 3> a/b, the optimal durations r^"^* and Tj''* can be explicitly given, 
and they are in continuity with the case without correlation p — 0.5. 



C. Other distributions of phases duration 

The model presented in section ITVl is minimal in the sense that the searcher has no memory. As seen in the previous 
section, a possibility is to add orientational memory. Another possibility is to add temporal memory. In the generic 
model, we have considered a "Markovian" searcher, in the sense that the rate of switching from one phase to the other 
is constant. It leads to an exponential distribution of the durati ons of the phases. In the following, we study the 
influence of the distribution of the duration of the phases (see also ( Chechkin et al.fa)l . l2009l : lMeerschaert et al. 
iTeiedor and Metzler 1 . 12009| ) for other types of correlations) . A first possibility is to study the effect of distributions 
which are peaked around the mean duration, or even d eterministic ("Bcnich ou et al.l [20071 ). A second possibility is to 
study the case of Levy-distributed "blind" phases as in iLomholt et al.i (j200^! 
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Figure 27 t°^^ as a function of 2(1 — p). Theoretical value of r{'*'* (|58p (dashed line) and theoretical value of Tj''* (|59p (solid 
lines), compared to the numerical minimization of the full exact mean search time, leading to tJ'^* (o) and Tj''' (□). a = 0.01, 
6=1 (green), a = 0.01, h = 100 (blue), a = 1, 6 = 100 (violet), a = 1, fe = 10"* (red), a = 100, 6 = 10"* (orange), a = 100, 
fe= 10*^ (brown), fc = 1, V = 1. 
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Figure 28 Comparison between the search without temporal memory (o) and the search with temporal memory (□). Static 
mode in two dimensions, tm. as a function of r2. k = 1, V = 1, b = 113, a = 10, ti = 2.6. 



1. Deterministic durations of the phases 

In the generic model previously described, we have considered exponential durations of phases, which correspond to 
searchers with no temporal memory. In the opposite case, a searcher with full memory skills could for example switch 
from one phase to the other at deterministic times instead of exponentially distributed times. The corresponding 
problem is not Markovian any more, which makes its analytical resolution much more complex. We present here a 
summary of a nume rical study of the effec t of such tempor al memory for a searche r with the static mode of detection 
in dimension 2 (see iBenichou et ah for details and iBenichou et al.l (|2006a[ ) for a semi-analytical treatment in 

1 dimension). First, this study shows that the optimal t°^^ and Tj^* are larger than in the case without memory, but 
are of the same order of magnitude (see figure [28l) . Second, such temporal memory decreases the mean search time. 
Indeed, a deterministic duration of the relocation phase avoids both the very short and very long relocations, that are 
inefficient. Third, and importantly, the gain from this temporal memory is quite low (less than 40% in an extended 
range of parameters, and decreasing with b/a increasing) as compared to the case with no memory (see figure [ 
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2. Levy distribution of the fast phase durations 

iLomholt et al.l ( 2008[ ) study analytically and numerically a 1-dimensional intermittent random walk whose duration 
of relocation phases is taken from a Levy law {p{l) oc with 1 < a < 2). Apart from this dis tribution of the 

duratio n of ballistic phases, this model is identical to the generic model presented in section lTvl fsee also lBenichou et all 
in the case of a diffusive mode of detection in one dimension (in the particular case of a point-like target 

The mean search time is evaluated with the exact formula (equation (9) of iLomholt et al 



2(ti +r2) 



OO 

DT.kl + 1 - A(fc„) 



(60) 



with kn = 2it/L, where L is the distance between targets , and A is the characteristic function of the distribution : 
X{k) = exp(— (T"|fc|") {p{k) = e^^'"^ p{x)dx) . The relation between cr, a, the velocity V and T2 (the mean duration 
of phase 2, which is defined since a > 1) is given by (equation (10) of ILomholt et al. j (|2008D 1 : 



2r(l - l/a)' 



(61) 



A mo re tractable approximate expression of the mean search time can be derived (equation (14) of ILomholt et all 
(I2008D 1 : 



{t) =2(Ti+r2) 



L 



4VSn 



V27rcr J 



C(«) 



(62) 



with C(a) = Yl^=i is the Riemann C function. 

As compared with the generic model of section HVl this model in troduces an extra par ameter a, which as could 
be expected enables a further minimization of the mean search time. ILomholt et al.l ( 20081 ) claim that Levy laws are 
more efficient than exponential laws because they have no second moment and therefore are not bound to the central 
limit theorem. 



Figure 29 tm as a function n, a being at the theoretical minimum (numerical minimization of equation [62]). lines : approached 
analytical formula [62l x : simulations without cut-off. + : simulations with cut-off at L. D = 1, V = 1. Left L — 10**, sienna : 
a = 1.6; red : a — 1.5; violet : a = 1.4; blue : a — 1.3; green : a — 1.2. Right L — 10^, sienna : a — 1.6; red : a — 1.5; violet : 
a = 1.4; blue : a — 1.3; green : a = 1.2. 



However, the Levy distribution is not the optimal distribution of the duration of ba llistic phases, and dist ributions 
with a finite second moment can perform even better, as opposed to what is claimed in ILomholt et al.l ( 20081 ). Indeed, 
relocations larger than the distance between two targets L cannot be profitable, and power law distributions are 
therefore inefficient in the regime of long times t > L/V. A simple example is given by a Levy distribution with 
an upper cut-off at L (see figure [I^ . For L ~ 10^, numerical simulations show that the optimum without cut-off 
is realized for approximately a ~ 1.4, with <,„ ~ 195 000: and the optimum with a cut-off at L is realized for 
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a ~ 1.3, with tm — 188 000, that is ~ 3.7% lower. For L = 10^, the optimum without cut-off is realized for 
a ~ 1.3, ~ 3 260 000; the optimum with a cut-off of L is reahzed for a ~ 1.2, — 3 060 000, that is ~ 6.5% 
lower. Truncated distributions with well defined second moment therefore outperform Levy distributions. Hence, 
intermittent random walks with Levy-distributed relocations decrease the mean search time more efficiently than in 
the case of exponentially distributed relocations, but it is not because of their infinite variance, and other distributions 

with second moment can perf orm even better . 

Another point discussed bv lLomholt et al" j (120081 ) is the robustness of the strategy : if L is misevaluated, the effi- 
ciency of the intermittent search with exponential relocation durations decreases more than for the Levy distribution. 
Truncated Levy is probably of intermediate robustness. 



D. The point of view of the target : Pascal principle 



In this review, we have addressed the question of determining optimal search strategies. One could also consider 
the opposite point of view and try to determine optimal survival strategies of targets. In the case where the target's 
motion is independent of the s earch er's motion, the response is actually given very simply by the so called Pascal 
principle (|Moreau et all l2003l |2004| ) for a broad class of situations . More precisely, suppose that the motion of 



the searcher is time and space homogeneous and that it satisfies the following property: starting from a position x 
(different from the target position), the transition probability to be at position y at time i > is always maximum for 
y — X. Assume that the target can perform any stochastic motion, independently of the searcher which is assumed 
to perform a stochastic motion. Then Pascal principle states that the survival probability of the target is maximum 
if the target remains immobile at its initial position. Of course, the validity of Pascal principle is restricted to special 
motions of the searcher: it holds for a diffusive motion, but not, for instance, if the searcher undergoes a ballistic 
motion with constant velocity. However, the validity conditions are satisfied if the searcher undergoes a ballistic 
motions with symmetrically distributed stochastic velocities, or if the displacements consist in teleportations, which 
are distributed symmetrically with respect to the initial position. In these cases, the best strategy for the target is to 
remain immobile. 



E. Other models of intermittent search 

In the past few years, several models relying on the mechanism of intermittent search have been developed in 
different contexts. We here brieffy review these models which are in essence similar to the generic case discussed in 
section IIVI and which broadens the field of application of intermittent search. 

a. \Oshanm etall \200± \200i ) 



2 2 




L L 



Figure 30 Model used bv lOshanin et al" j (|2009l . l2007l ). 



lOshanin et al. ] Hooilloo^) propose a model very similar to the diffusive mode in one dimension of the generic 



model, but in discrete space, on a infinite lattice. At each time step, with probability a, the searcher jumps to the 
neighboring node of the line (with equal probabilities for each side, which corresponds to diffusion). With probability 
1 — a, it stays off-lattice during a time T and after this time, it lands at a distance L from its initial position (once 
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again, with equal probabilities for each side) (see figure 150)1 . This phase is equivalent to a ballistic non- reactive phase. 
Its duration is exactly T, whereas the duration of the diffusive phase with target detection is exponentially distributed, 
with mean duration 1/(1 — a). There is one target, but an infinite set of searchers, initially randomly distributed. 
The quantity maximized is the probability that at a given time t, the target has already been found by any of the 
searchers. Oshanin et al. do find an optimal a, but dependent on t. If L and T are fixed, then the optimization 
with respect to a leads to a°^* oc t^^^, which can be very small. If T only is fixed , the optimization with respect 
to both a and L gives : a°P* = 0.5 and L°p* ^ \/t/\n{t). Note that in the case T = 1, at the optimum the time 
has to be equally shared between the two phases, which is reminiscent of the result obtained in the framework of the 
simple model of facilitated diffusion (section IIII.A.2p . Last, if V ^ L/T is fixed, it is found that L°p^ ^ \/t/\n{t), 
a°P* ^ 1 — AV \n{t) / {3\/t) , which means than the mean duration of the nearest-neighbors ph ase is 3 /4 the du r ation o f 
a large move. This model can actually be seen as the lattice version of the model given in ( Benichou et al.l . l2006af l. 



where similar results were obtained : (i) existence of a global minimum of the search time ; (ii) the ratio of times 
spe nt in both phases at th e minimum is given by a numerical constant ; (iii) this numerical constant is equal to 1/2 
in ( Benichou et al.l , [2006a[ ) (where the time spent off the lattice is exponentially distributed) instead of the 3/4 given 



above (for deterministic times off lattice) . 
b. \Roio et aJ.I f200£l ) 




target 




Figure 31 Model used bv lRoio et ail (120091 ). 



iRojo et all (|2009t ) propose a model which displays some similarities with the previous model (see figure [31]) . The 



search domain is also a one-dimensional discrete infinite lattice with one target, there are also an infinite set of 
searchers, and the quantity optimized is also the probability that the target is found by any of the searchers at a given 
t. The detection phase consists of jumps to the nearest-neighbors, with a given frequency. Such a rule is equivalent 
to diffusion. The non-reactive phase consists of jumps to the next nearest neighbors. It is again diffusion, but if the 
jump frequency is the same as in the other phase, it is a faster diffusion. In both phases, there is a fixed rate of 
switch to the other phase, leading to exponentially distributed durations of the phases. If one of the mean durations 
is fixed, the probability that the target is already found at t is minimized for a finite duration of the other phase. But 
the optimum is for infinitely short phases, enabling the searcher to combine the faster diffusion of one phase and the 
detection capacities of the other phase. 

c. iReingruber and Holcmail t200(\ 



iReingruber and Holcmann 20091 ) propose a model which is also diffusive/ diffusive (see figure |32|) . They study this 



model first in one dimension : the searcher's starting point is at one extremity of a segment, a reflecting boundary. 
The target is at the other end of the segment. However, in phase 1 (diffusion of coefficient Di), the target can 
be found, whereas in phase 2 (diffusion of coefficient D2), both extremities are reflecting. There are fixed rates of 
switching from one phase to another. The results show that there are two regimes : ii Di > D2, straightforwardly, 
the optimum for the searcher is to be in phase 1 only; if D2 > Di, the optimum is to switch very rapidly between 
the two phases, such as to spend almost all the time in the faster phase 2, but not to miss the target. This model 
is extended to a 3-dimensional ball (the initial position is almost without importance in this geometry), but with 
a target of radius a on the border ( which is reflecting everywhere else). The two phases are defined like in one 
dimension. IReingruber and Holcmanl (|2009t) give two limits in this case. I t can be noticed that t he expression we have 
obtained in the generic model for the diffusive mode in 3 dimensions (see iLoverdo et al. j lloogB)) could be used, with 



D2T2. In fact, our calculations use a "diffusive/diffusive" approximation, with an effective D2 = 3V T2. The 
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optimization will be quite different, because the dependence in T2 is dramatically changed if instead of a fixed 1)2, D2 
is a function of T2. Indeed, the optimum for our generic model is for finite t\ and T2, whereas, even if not explicitly 
calculated, it is probable that the optimum for diffusion/diffusion in three dimensions is similar to the one-dimensional 
case, i.e. for phases durations as small as possible. The goal of this model is to study cellular signaling, with a ligand 
binding to a target which will transmit a signal. 



d. iBressJoff and JVewb-d if200£l ); lJVewbv and Bressto^ if200£l) 

V+ V- 



target 
I 




Figure 33 Model used by IBressloff and Newbvl (|2009l '): iNewbv and Bresslofj l|2009r ). 



IBressloff and Newbvl ( 20091 ) present another model applied to intracellular transport, more precisely here to the 



transport of mRNA granules inside neurons. They present a model in one dimension, standing for example for an 
axon with little branching. The starting point is at one extremity, which is reflecting : it models granules produced in 
the soma of the neuron and that have to be exported to the axon. The target, a synapse, is somewhere in the segment. 
The other end of the segment is an absorbing boundary, representing that the vesicles containing the mRNA can be 
degraded, or that there can be other targets further away in the axon that can absorb the searcher. To complete 
this idea that there are several targets that are not equivalent, and that these targets are in competition, they also 
calculate explicitly the probability that the searcher finds a target more often than the others. In this model, there 
are 3 states (see figure [33|) : an immobile detection phase, similar to the static mode, switching to ballistic modes 
with probability a per unit time; a ballistic phase in direction +, with speed and with a transition rate to the 
detection mode a ballistic phase in direction — , with speed and with a transition rate to the detection mode 
/3_. During the two ballistic phases, the searcher cannot detect the target. Movement is biased to the direction + if 
>z;_//3_. 

The results are based on the fact that on the segment, there are two contradicting constraints : maximizing the 
hitting probability (as the searcher can be degraded before finding the target), and minimizing the time to find the 
target when the target is found. Indeed, if there is more bias, the target will be missed more often, but when found. 
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the search time wiU be smaUer. With a fixed hitting probabiUty, the mean first passage time to the target (on the 
condition that the target is found) is minimized when there is more bias. In other words, unidirectional motion is 
b etter than bi-directional mo tion in this case. 

iNewbv and Bresslofil ( 20091 ) extend this problem to the case of a directed tree. In this case, unidirectional motion has 
a drawback : a wrong branch can be taken, annihilating any possibility to find the target. Biased bidirectional motion 
can be seen as an effective combination of a ballistic and a diffusive motion. It exists a critical hitting probability p* . 
If the mean first passage time to the target is minimized given that the probability of finding the target is a given 
p < p* , unidirectional motion is better; but if the given probability is p > p*, there is an optimal finite bias which 
minimizes the mean search time in case of success. 



e. \RamezanDoui if2007l) 




log(k) 



Figure 34 Model presented bv iRamezanpouJ (j2007l ). 



Intermittence in networks such as the tree we have just seen is an interesting extension. iRamezanpouil (|2007t ) 
proposes (see figure 15^ to e xplore a network in which t he degree (= number of neighbors) distribution is p{k) cx k~^, 
constructed as proposed by iBarabasi and AlbertI (|l999l ). or with some modifications. On this finite network, at each 
time step, the searcher chooses randomly one of the edges connected to the node where it is, and goes to the node 
connected by this edge. Every tw, the searcher jumps to a completely random node. The question is whether the 
mean time to cover the nodes and the edges of the network can be optimized as a function of t^,. For the nodes, the 
random jumping is a way to visit all the nodes with equal probability, thus should be as small as possible. For the 
edges, there is an optimal finite i^,. Indeed, if is small, most edges visited will emanate from low-connected nodes 
(as the low connected nodes are the more numerous nodes, such edges are more likely to be visited after a random 
jump), but if is large, the searcher would spend most of its time on the edges connecting high degree nodes, and 
will take time to explore the whole network, especially for remote edges connecting nodes of low degree. 



F. Designing efficient searches 

As seen previously, intermittent reaction paths are involved in various search problems involving biomolecules at 
the microscopic scale, as well as biological organisms at the macroscopic scale. Simple analytical models show that 
intermittent transport can actually minimize the search time. A reason why such intermittent trajectories are widely 
observed could be simply that they constitute very generic optimal search strategies, and consequently they could 
have been selected by evolution. 

Beyond modeling what is observed in real-life biological examples, such intermittent strategies could also be used 
to design searches, at the microscopic and macroscopic scales. We briefly discuss here potential applications at the 
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microscopic scale (for more details see iBenichou et aL j (l2008ari 'l 






Without flow 



With a flow 



Figure 35 Design of heterogeneous chemical reactions, with targets (green disks) here fixed on a 2-dimensional surface. The 
reactant either diffuses in the volume (blue), or diffuses on the surface (orange thick line). The flow is represented by magenta 
arrows. 



Heterogeneous chemical reactions, where the reactive targets are located at an interface, either 1-dimensional 
(polymer) or 2-dimensional (surface) are intrinsically intermittent. Indeed, the reactants can either diffuse in the bulk 
volume, where target cannot be found, or bind to the interface and diffuse more slowly (see figure [35l left). 

Beyond obvious optimizations (increasing the target and the reactant concentrations, increasing the diffusion coef- 
ficients of the reactant in the bulk or at the interface, etc.), the mean durations of the phases (free or bound to the 
interface) are the main adjustable parameters enabling to minimize the mean search time, and therefore to maximize 
the reaction rate. The main idea is that in the "teleportation" approximation (see section |III.A.2[) . after a bulk 
excursion, the distance between the reactant landing point on the surface and its starting point is larger than the 
typical distance between targets. In such a regime, each new bound phase is independent from the previous one, and 
the trajectory overlap is limited, which enhances the reactivity. 

The mean durations of the phases can actually be tuned in real systems. First, the mean time a reactant remains 
bound to the interface depends on its affinity with the interfac e, which could be tuned. Second, the mean time 
spent in the bulk is main ly controlled by the confinement volume ( Benichou et al.l , l2005al : iBlanco and Fournied . [2003l; 
ICondamin et al. ], ioozB). The confinement volume has therefore to be as small as possible, but it should be large 
enough to make the "teleportation" approximation valid. This constraint defines a critical volume, and can actually 
be by-passed by applying a hydrodynamic flow parallel to the surface, which makes the teleportation approximation 
valid even for a very short bulk excursions, provided th at the velocity of the flow is high enough (see figure 1551 right). 
In this regime it can be shown (jBenichou et al 

lioosi) 

that the reaction rate can be optimized by tuning the affinity 
of the reactant for the interface in a similar way as in section IIII.A.2I 

At the molecular level, we stress that intermittent transport could also be useful for in vitro chemistry. Indeed, 
we have shown that intermittent transport naturally pops up in the context of reaction at interfaces, where reactants 
combine surface diffusion phases and bulk excursions, and could permit to enhance reactivity. In this case, adjusting 
chemically the typical association time of the reactants with the interface, makes possible to optimize the reaction 
rate. 
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VI. CONCLUSION 

Intermittent search strategies rely on a simple mechanism : the searcher alternates between two phases, one during 
which the target can be detected, but with slow motion, and another of faster motion but without target detection. 

This mechanism of intermittence has emerged from the observation of real-life biological searches at various scales. 
At the macroscopic scale, an example is given by animals searching for hidden food, which alternate between fast 
ballistic relocation phases with no target detection, and phases of slower motion aimed at detecting the target. A 
simple model based on this observation permits to show analytically that the mean search time can be minimized as 
a function of the phases mean duration. There is one single way to share time between the two phases to find the 
target as fast as possible. This intermittent search is then an optimal search strategy. With this respect, this model 
is an alternative to the famous Levy walks model which is optimal only in restrictive conditions. 

Intermittence is also observed at the microscopic scale. Indeed, for some biochemical reactions in cells, which involve 
a very low concentration of reactants, reaction pathways are not always simple Brownian trajectories. They can rather 
be qualified as intermittent, since they combine slow diffusion phases on the one hand, and a second mode of faster 
transport on the other hand, which can be either a faster diffusion mode as in the case of DNA-binding proteins, or 
a ballistic mode powered by molecular motors in the case of intracellular transport. Analytical models actually show 
that such intermittent trajectories are very efficient, since they significantly reduce reaction times. Interestingly, it 
is shown that reaction rates can even be maximized by adjusting simple biochemical parameters. The gain is small 
in dimension three, but for lower dimensional structures such as membranes (2D) or polymers or tubular structures 
(ID), the gain can be very large at low target concentration. Such efficiency - and optimality - could explain why 
intermittent transport is observed in various forms in the context of reactions in cells. 

Since these intermittent search strategies are observed at various scales, one could suggest that they constitute a 
generically efficient search mechanism. The systematic analysis of a generic model in the framework of intermittent 
random walks, in 1, 2 and 3 dimensions, and for three different descriptions of the slow reactive phase, permits to 
assess quantitatively the robustness of this mechanism. In fact, this study shows that the optimality of these search 
strategies is a widely robust result. Finally, if intermittent random walks are observed in real biological systems at 
various scales, it is probably because they do constitute an efficient search strategy. Beyond these modeling aspects, 
one can suggest that such intermittent strategies could also be used to design optimizable search strategies. 
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Appendix A: Recap on random walks and Levy processes 

Regular random walks obey Gaussian statistics, and have a mean square displacement growing linearly with time 
■ (''^(^)) ~ ^"i with a ~ 1. Inversely, transport processes characterized by non- linear scalings with time of the mean 
square displacement are termed "anomalous", either subdiffusive if a < 1 or superdiffusive if a > 1. In this review, 
we shall make use of standard models of subdifFusion (Continuous Time Random Walks, Diffusion on Fractals and 
Fractional Brownian Motion) and supe rdiffusion (Levy flights and Lev y walks), whose definitions are briefly reminded 
here for consistence (see for instance iBen-Avraham and HavlinI (|2000() for a more complete discussion) . 



1. SubdifFusion 

a. Continuous time random walks 



A first class of models leading to subdifFusion stems from continuous time random walks (CTRWs) and their 
continuous space limit described by fractional diffusion equations. The anomalous behavior in these models originates 
from a heavy-tailed distribution of waiting times : at each step the walker lands on a trap, where it can be trapped 
for extended periods of time. Technically, the CTRW is a standard random walk with random waiting times, drawn 
from a probability density function 'ip{t). The CTRW model has a normal diffusive behavior if the mean waiting time 
is finite. For heavy-tailed distributions such that 

V'(i) oc -j^, at large times, (Al) 

the mean waiting time diverges for /3 < I and the walk is subdiffusive with a = j3. 

When dealing with a tracer particle, traps can be out-of-equilibrium chemical binding configurations, and the 
waiting times are then the dissociation times; traps can also be realized by the free cages arou nd the tracer in a dense 
hard sph ere-like crowded environment, and the waiting times are the life times of the cages (jCondamin et al 1 120081: 
ISaxtonl . [l99& . ,2007h . 



b. Diffusion on fractals 



Another kind of model for subdifFusion relies on spatial inhomogeneities as exemplified by diffusion in deterministic 
(such as Sierpinski gasket) or random fractals (such as critical percolation clusters). The subdiffusive behavior is in 
this case caus ed by the presence of fixed obstac les that create numerous dead ends, as illustrated by De Gennes's ant 
m a labyrinth (|de Gennes l . ll976l:ISaxton I . ll994[) . This results in an e ffective subdifFu s ion in the embedding space, with 
an ex ponent a < 1, who se value depends on the fractal structure ( d'Auriac et al.l , Il983l : iBen-Avraham and HavlinI . 
120001 : (Bunde and HavlinI . 1 19911) . 



Fractional Brownian motion 



Fractional Brownian motion (FBM) is a third model of subdiffusion ( Mandelbrot and van Nessl . Il968l) . usually 



defined for systems in dimension 1. It was introduced to take into account correlations in a random walk: the state of 
the system at time t is influenced by the state at time t' < t. More precisely, it is a Gaussian process of autocorrelation 
function of the form 

{X{h)X{t2)) oc t\" + - \h - 12^, (A2) 

with < iJ < 1/2, so that a = 2H < 1 (FBM can also be defined for 1/2 < iJ < 1, but in this case it leads to 
superdiffusion). Note that Brownian diffusion is recovered for H = 1/2. FBM is used to describe the motion of a 
monomer in a polymer c hain or single file diffusion. Recently, it has also been proposed to underline the diffusion in 
a crowded environment ( Szvmanski and Weissl |2009|) . 
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2. SuperdifTusion 

a. Levy flights 

Levy flights are random walks such that, at each step t, the walker jumps in some random uniformly distributed 
direction, to a distance r drawn from a probability density function 

P(^) °^ (A3) 
It can be shown that, if /3 < 2, superdiff'usion emerges, with a = 2//3, while, if /3 > 2, regular diffusion is recovered. 



b. Levy walks 



Levy walks differ fro m Levy flights in that, now. the time to make a step of size r is taken proportional to r 
(jShlesinger et al~ I [l983). Physically, it can be seen as a random walker performing jumps still drawn from the 
probability density function 



p{r) (X 



,1+13 ■ 



(A4) 



but this time at a constant velocity. The resulting mean square displacement can be shown to be given by 





if < /3 < 1, 




if/3 = l. 




if 1 < /3 < 2, 


tint 


if /3 = 2, 




if /3 > 2. 



(A5) 



Applications of Levy walks are given in the main text, in the context of random search problems. Note that, in 
this context, the terms "Levy walks" and "Levy flight" are often used indifferently to design "Levy walks". 

Appendix A: Mean first-passage times of intermittent random walks 
1. Dimension one 



a. Static mode 



In this section we assume that the detection phase is modeled by the static mode. Hence the searcher does not 
move during the reactive phase 1, and has a fixed reaction rate k per unit time with the target if it lies within its 
detection radius a (see figure [21]). It is the limit of a very slow searcher in the reactive phase. 



a. Equations 

Outside the target (for x > a), we have the following backward equations for the mean first-passage time : 



V-^ + -{ti-t^ 

ax T2 



= —1, and 



(A2) 



54 



Inside the target (x < a), the first two equations are identical, but the third one is written 

1 it + f7 /I 



n 2 



+ fc tl = -1. 



We introduce t2 ~ *^ and = *^ . Then outside the target we have the fohowing equations 

ax T2 

1 



(t2-ii) = -l. 



Inside the target the first two equations are identical, but the last one writes : 

—t2 -(— + k]ti = -1. 



(A4) 

(A5) 
(A6) 
(A7) 

(A8) 



Due to the symmetry x —x, we can restrict the study to the part x G [0, a] and the part x € [a, b]. This symmetry 
also implies : 



dx 



0, 



dx 



x=0 



x—b 



In addition, continuity at x = a for and gives: 



j_d,in / \ id.out / \ 

t2 — aj ~ C2 — '^j- 



(A9) 
(AlO) 

(All) 
(A12) 



This set of linear equations enables us to explicitly determine ti, t2, 4 inside and outside the target. 



b. Results 

An exact analytical expression of the mean first passage time to the target is then given by : 

(6-a)3 /3(6-a)2 



T1+T2 ( b 



fcri 3V^V| Vt2 



-coth 



VT2I3 



{b - a)l3' 



where P = ^(fcri)-! + 1. 
We obtain in the limit b':$> a 



tm = in + T2) 



f p fj_ \ b 



We use this approximation (jA14p to find ti and T2 values minimizing t„i : 



opt 



alb 
Vk \ mi 



1/4 



(A13) 



(AM) 



(A15) 



opt _ / b 



(A16) 



Importantly, the optimal duration of the relocation phase does not depend on fc, i.e. on the description of the detection 
phase. 
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1 / Opt 5 



Figure 36 Static mode in one dimension. Exact expression of trr. CATsIi (1 ines) compared to tlie approximation of tm (|A14[) 
(symbols), both rescaled by C*- from (fUS)) . t^"- from (fXTHl) . 1/ = 1, fc = 1. 6/a = 10 (green, □), fe/a = 100 (red, o), 
6/a= 1000 (blue, +). 



b. Diffusive mode 



We now turn to the diffusive modeling of the detection phase. The detection phase 1 is now diffusive, with immediate 
detection of the target if it is within a radius a from the searcher (see figure I^Tj) . 



a. Equations 

Along the same lines, the backward equations for the mean first-passage time read outside the target {x > a) 



,dtt 1 



ax T2 



(A17) 



dx 



+ -(ti-i^) = -l, 

T2 



dHi 



(A18) 
(A19) 



and inside the target [x < a) 



dx 



T2 



(A20) 
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- -t- = -1, and (A21) 

ax T2 



ti = 0. (A22) 
Boundary conditions result as previously from continuity and symmetry. 

b. Results 

Standard but lengthy c alculations lead to an exact expression of the mean first detection time of the target tm given 
in ( Loverdo et al Three regimes can be identified. 



In the first regime (6 < ■^) intermittence is not favorable. 



• For b > y and -^iy^ < 1 intermittence is favorable. In the low target density (6 ^ a) limit, the following 
approximation of the mean first passage time around its minimum can be obtained : 



This expression gives a good approximation of tm in this regime, in particular around the optimum ((see 
figure |37)) ) . The simplified tm expression (jA23|) is minimized for : 



opt _ sj 2b^D 

'^2 - \/ gyi ' 



(A25) 



24 DV^ 

This compares to the case without intermittence according to 

opt tdis _ ,f¥ fbV\^ ^.^fbV 



C'-a/St^- (A26) 



• For b > y and 1 ^ '^iyj intermittence is favorable. 
In this regime, the mean search time can be approximated by: 

b / a ab 



tm ^ -(ri + T2) ( _ , ^ + :;775-2 ) ■ (A28) 



This expression gives a good approximation of tm-, at least around the optimum (see figure I38p . which is 
characterized by: 



opt _ Db 



48V^a 



(A29) 



rr - f (A30) 

I aV pb 

nam ~ — = \ —■ (A32) 

^ 2^/3 V a ^ ' 

Note that the gain can be very large at low target density, and that the value obtained for Tj''* is the same as 
in the static mode. 
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Figure 37 DiflFusive mode in 1 dimension, -f^, exact expression (line), and approximation in the regime of favorable inter- 
mittence and > 1 (|A23|) (symbols). a = l and 6 = 100 (green, o), a = 1, 6 = lO" (red, +), a = 10, b = 10^ (blue, □). 

D = 1, V = 1. Tj''* is from expression (|A24p . Tj''* is obtained from expression (|A25|) . 




Figure 38 Diffusive mode in 1 dimension. -[^, exact expression (line), and approximation in the regime of favorable inter- 

mittence and < 1 (fA28)l (symbols), a = 10 and 6 = 100 (green, o), a = 10, 6 = 1000 (red, +), a = 100, b = lO" (blue, 

O). D — 1, V = 1. T°^* is from expression (fl29)) . T°^^ is obtained from expression ()A30p . 



c. Ballistic mode 



We now treat the case where the detection phase 1 is modeled by the ballistic mode (see figure [2T|) . This model 
schematically accounts for the general observation that speed often degrades perception abilities. Our model corre- 
sponds to the extreme case where only two modes are available : either the motion is slow and the target can be 
found, or the motion is fa st and the target cannot be found. Note that this model can be compared to the model of 
IViswanathan et aP ( 1999[ l. where there is only the detection phase. 
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a. Equations 

The backward equations read outside the target (x > a) 



+ + (A34) 



dx n V 2 2 



T2 \ 2 

Inside the target (a; < a), one has t^'™(x) = t^'"^{x) ~ 0, and : 

- -*2"'" = -1- (A38) 

b. Results 

In the case where the phase 1 is modeled by the baUistic mode in one dimension, we have calculated the exact mean 
first passage time i„j at the target. can be minimized as a function of ti and T2, yielding two possible optimal 
strategies : 

• for vi > vf — V^y^, intermittence is not favorable : t°^* — >■ 00, t^^* — > 0, 

• for vi < vf — V^^^, intermittence is favorable, with r°^* — > and t^'''^ sa 'W\f^ ■ "^^^ gain reads 

gam « — W-. (A39) 

2vi\b 

This shows that the gain is larger than 1 for vi < vf — V^y^^, which defines the regime where intermittence 



^'i - " —V b 
is favorable. 



Note that the model studied bv IViswanathan et al. shows that when targets are not revisitable, the optimal 



strategy for a single state searcher is to perform a straight ballistic motion. This strategy corresponds to ti — > 00 
in our model. Our results show that if a faster phase without detection is allowed, this straight line strategy can be 
outperformed. 

d. Conclusion in one dimension 

Intermittent search strategies in one dimension share similar features for the static, diffusive and ballistic detection 
modes. In particular, all modes show regimes where intermittence is favorable and lead to a minimization of the search 
time. Strikingly, the optimal duration of the non-reactive relocation phase 2 is quite independent of the modeling of 

the reactive phase : T2 ^* = ^ for the static mode, for the ballistic mode (in the regime vi < vf o:^ ^ \f^) ' 
for the diffusive mode (in the regime b > y and F v a)- This shows the robustness of the optimal value T2^^ . 
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2. Dimension two 

a. Static mode 

We study here the case where the detection phase is modeled by the static mode : the searcher does not move 
during the detection phase and has a finite reaction rate with the target if it is within its detection radius a (see 
figure m]). 



a. Equations and results 



ft 




Figure 39 Static mode in two dimensions. Simulations (symbols) and analytical approximate (|A43[l (lines), k — 1, V — 1, 
& = 56; a = 10 (red, +) (r"*" = 2.41, r""' = 11.2), a = 1 (blue, o) (t""' = 0.969, r^* = 1.88), a = 0.1 (green, □) (r""* = 0.348, 
Tj''' — 0.242). Left : mean search time tm as a function of r2/a, with ri — t°^^ . Right : mean search time tm as a function of 

n, with T2 = 



The mean first passage time at a target satisfies the following backward equations ( Redneil . l200l[ ): 

— rihii^) - h{i^))de^ - kia{-^)t^{i^) = -1. (A40) 



2t:t- 



1 Jo 



f ■ Vri2(^) - -iHl^) - h{l^)) = -I (A41) 

The function is defined by laC?^) = 1 inside the target (if |~?^| < a) and la(^) = outside the target (if \'^\ > a). 
In the present form, these integro-difFerential equations (completed with boundary conditions) do not seem to allow 
for an exact resolution with standard methods. t2 is the mean first passage time to the target, starting from 1^ in 
phase 2, with speed 1^, of angle 8-v, and with projections on the axes Vx, Vy. i and j can take either x or y as a 
value. The following decoupling assumption is introduced : 

{V,Vjt2)e^ ^ {V^V,)e^{t2)e^ (A42) 

and leads to the following approximation of the mean search time, which can be checked by numerical simulations 
(see figure [39]) : 

. + kr^Xy^ - x^r^^^ + \ [Sy^ + (1 + fcn) (V HvM + {y^ - ^')(^' - + 8))] } , (A43) 



where x = a/— ^^r- — and ii = a/— ^^r- — ■ (A44) 

V 1 + A:ti Vt2 V 1 + fcn Vt2 ^ ' 

In that case, intermittence is trivially necessary to find the target. In the regime b ^ a, the optimization of the search 
time (|A43|) leads to : 



opt 

^1 -\VkJ 



/_t^\i/2 / 21n(Va)- l y/^ 



(A45) 
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rr = ^{Hh/a)-l/2f'\ (A46) 
and the minimum search time is given in the large volume limit by : 

opt_^_ 2^/^ (g^ - 4&^) hi(Va) + 26^ - 
^""a^fc Vyfc^ (21n(6/a)- 1)3/4 

72 (960^62 _ 192^4) in2(5/a) + (19264 _ 144^2^2) ij^(^/^)_^4g^2^2_ 47^4 ^4 
~ 48a62y (21n(fe/a)- 1)3/2 ■ 



b. Diffusive mode 



We now assume that the searcher diffuses during the detection phase (sec figure [2T| . For this process, the mean 
first passage time to the target satisfies the following backward equation (jRedner, . 



27rTi Jo 



(A48) 



T2 



(A49) 



with ii(~P") — inside the target (r < a). The same decoupling assumption as for the static case sis used (jA42|) . 
It eventually leads to the following approximation of the mean search time, which can be checked by numerical 
simulations (see figure |40| : 



t,n (n + T2) ^ — -2 { aa(b /a 



with : 



{ niii M L_ a^DTll■i-A\T^{b a)W-Aa'h'+a''^ 
I aa{V - I)— ^ — ^ '' ^ \ , 



(A50) 



with L± —Iq 



(Ii(6a)Ki(aa) - Ii(aa)Ki(6a)) 



± a\l Dt2 Ii 



(Ii(6a)Ko(aa) + Io(aa)Ki(fea)) 



(A5I) 



and : 



M =1, 







(Ii(6a)Ko(aa) + Io(aa)Ki(&a)) 



^Dt2, 

^ aV^^2 

a[b'^ — a?)'- 



(Ii(&a)Ki(aa) - Ii(aa)Ki(6a)) , 



(A52) 



where a — {\/{Dti) + 1/{Dt2)Y^'^ and D = v^T2. The minimization as a function of ti and T2 is as follows. 



a. a <b <^ D/V : intermittence is not favorable 

In that regime, intermittence is not favorable. Indeed, the typical time required to explore the whole domain of 
radius b is of order b^/D with diffusive motion, which is shorter than the corresponding time b/V with ballistic motion. 
As a consequence, it is never useful to interrupt the diffusive phases by mere relocating ballistic phases. T he mean 
first passage time to the target in this optimal regime of diffusion only is obtained using standard methods ( Redneij . 
1200 ll) and reads in the limit b ^ a : 

tdiff = (-3 + 4 hi -) . (A53) 
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Figure 40 Diffusive mode in two dimensions. Simulations (symbols) versus analytical approximate (jA50|) (line) of the search 
time, rescaled by the value in the absence of intermittence tdiff as a function of T2 (left) and In(ri) = (right), for D = 1, V = 1, 
b = 226. Left : a = 10, n = 1.37 (green □); a = 1, n = 33.6 (blue o); a = 0.1, n = 213 (red +). Right: a = 10, T2 = 15.9 
(green, □); a = 1, r2 = 13.7 (blue, o); a = 0.1, T2 = 22 (red, +). 



b. a <^ D/V <^ b : £rst regime of intermittence 

In this second regime, one can use the following approximate formula for the search time: 



'™ - ADV^a^ VJ' {'^"^^/"^ ^ ' 2i^(ln(aa) + 7 - ln2) ^ (A54) 



7 being the Euler constant. An approximate criterion to determine if intermittence is useful can be obtained by 
expanding t„i in powers of 1/ti when ti — )■ cxj (ri — > 00 corresponds to the absence of intermittence), and requiring 
that the coefficient of the term 1/ri is negative for all values of T2 . Using this criterion, we find that intermittence is 
useful if 

exp(-7/4 + -i)Vb/D - 4 ln(6/a) + 3 > 0, (A55) 
In this regime, using Ea (jA54[) . the optimization of the search time leads to : 



opt 4 In w — 5 + c opt b y/4\iiw — 5 + c 



where w is the solution of the implicit equation w = 2Vbf(w)/D with : 

-8(lnuj)2 + (6 + 81n(6/a)) \nw- 101n(5/a) + 11 - c(c/2 + 21n(a/6) - 3/2), (A57) 



V41nw — 5 + c 2 



and c — 4(7 — ln(2)), 7 being the Euler constant. An useful approximation for w is given by : 



D ' \2Dhi{b/a) 

The gain for this optimal strategy reads : 

tdiff 1 41n6/a- 3 + 4aV6^ - aV&* f 1 wD 41nw-7 , , , 



C* 241nVa-3 + 2(41nw)ln(6/au;) V41nw-5 (4 In u> - 5)3/2 

If intermittence significantly speeds up the search in this regime (typically by a factor 2) , it does not change the order 
of magnitude of the search time. 
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Figure 41 Ballistic mode in two dimensions. ln{tm) as a function of ln(r2). Simulations (symbols), diffusive/diffusive approxi- 
mation (|A50|) with HA65|) (colored lines), ri — limit (IA66|I (black line), n — ^ oo (no intermittence) (|A64|) (dotted black line). 
6 = 30, a = 1, V = 1. n = (black, ★), n = 0.17 (yellow, •), n = 0.92 (green, o), n = 5.0 (blue, x), n = 28 (purple, o), 
ri = 150 (red, +), ri = 820 (brown, □). 



c. D/V <^ a <^ b : "universal" regime of intermittence 

In the last regime D/V <C a <C &, the optimal strategy is obtained for : 



D ln^(&/a) a 

'1 — : 

and the gain reads : 



tdiff V2aV ( 1 lo (2/V21n(6/a)-l) 



41n(6/a)-3i^ (2/ ^2 In(Va) - l) 2V21n(Va) 



Here, the optimal strategy leads to a significant decrease of the search time which can be rendered arbitrarily smaller 
than the search time in absence of intermittence. 



c. Ballistic mode 



In this case, the searcher has access to two different speeds: one {V) is fast but prevents the searcher from finding 
its target, and the other one (uj) is slower but enables the searcher to detect the target (see figure [^T|) . 



a. Simulations 

Since an explicit expression of the mean search time is not available, a numerical study is performed. Exploring 
the parameter space numerically enables to identify the regimes where the mean search time is minimized. Then, for 
each regime, approximation schemes are developed to provide analytical expressions of the mean search time. The 
numerical results presented in figure [41] suggest two regimes defined according to a threshold value vf of vi to be 
determined later on : 

• for vi > vf^ t,„ is minimized for T2 — > 

• for vi < uf , tm is minimized for ti — > . 



b. Regime without intermittence (t2 — 0, ri — >■ 00 j 

Qualitatively, it is rather intuitive that for vi large enough (the precise threshold value vf will be determined next), 
phase 2 is inefficient since it does not allow for target detection. The optimal strategy is therefore r2 — )- in this case. 
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In this regime, the searcher perforins a ballistic motion, which is randomly reoriented with frequency 1/ri . Along 
the same line as in IViswanathan et alJ ( 1999[ l (where however the times between successive reorientations are Levy 



distributed) , it can be shown that the optimal strategy to find a target (which is assumed to disappear after the first 
encounter) is to minimize oversampling and therefore to perform a purely ballistic motion. In our case this means 
that in the regime T2 — > 0, the optimal ri is given by r^*^* — >■ 00. 

In this regime, we can propose an estimate of the optimal search time thai- The surface scanned during St is 2avi5t. 
p{t) is the proportion of the total area which has not yet been scanned at t. If we neglect correlations in the trajectory, 
p{t) is solution of : 



Then, given that p{t = 0) = 1, we obtain : 



$ = -^^. (A62) 
at TTO^ 



p(t)^cxp(-^L (A63) 



and the mean first passage time to the target in these conditions is : 

tbai = -~ / t-^dt=- . (A64) 

Jo at 2avi 

This expression yields a good agreement with numerical simulations. Note in particular that ti,ai oc — . 



c. Regime with intermittence ri — >■ 

In this regime where vi < vf, the numerical study shows that the search time is minimized for ri — >■ (see figure Hi]) . 
We here determine the optimal value of T2 in this regime. To proceed we approximate the problem by the case of a 
diffusive mode previously studied (jASOp . with an effective diffusion coefficient : 

D=^-tl. (A65) 

This approximation is very satisfactory in the regime ti (see figure HT|) . 

We can then use the results of the previous section for the diffusive mode in the ri — )■ regime and obtain: 

^ 1 (3 + 41n(f))6'i-4a2fe2 + „ ^ , v.-/. , ^ 

The calculation of Tj''* minimizing t,„ then gives: 





a 



b\ 1 



^r = ^JlM7)-:^- (A67) 



Finally the gain reads 



d. Determination of vf 

Note that an estimate of vf can be obtained from (jA68|) as the value of vi for which gam — 1 : 

.F.^flnf^^y"^^^. (A69) 

It is noteworthy that intermittence is less favorable with increasing b. This effect is similar to the 1 dimensional case, 
even though it is less important here. It can be understood as follows: at very large scales the intermittent trajectory 
is reoriented many times and therefore scales as diffusion, which is less favorable than the non intermittent ballistic 
motion. 
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d. Conclusion in dimension two 

Remarkably, for the three different modes of detection (static, diffusive and ballistic), we find a regime where 

intermittence minimizes the search time for one and the same Tj^*, given by Tj^* = -^y^ln (^) — i. As in one 
dimension, this indicates that optimal intermittent strategies are robust and widely independent of the details of the 
description of the detection mechanism. 



3. Dimension three 

a. Static mode 

We study in this section the case where the detection phase is modeled by the static mode, for which the searcher 
does not move during the detection phase and has a finite reaction rate with the target if it is within a detection 
radius a (see figure [^T|) . 



a. Equations 

Denoting ti{r) the mean first passage time to the target starting from a distance r from the target in phase 1 
(detection phase), and t2,e,4>{r) the mean first passage time to the target starting from a distance r from the target 
in phase 2 (relocation phase) with a ballistic motion in a direction characterized by 9 and 0, we obtain : 



V.\/t2,e,4, H (ti - t2,e,4>) = -1- 

T2 



Then outside the target [r > a) : 



and inside the target (r < a) 



Ti \ATr Jq 



— — / d6'sin6' 



J d9 sin 9 J d<pt2,e,ci, - hj = -1, 



— (— I d9sm 



Ti 47r J„ 



2tt 



d(j)t 



2,0,(1} 



— + k]ti 



Defining t2 ~ d9sm9 j^^^ d(j)t2,e,if>, one obtains outside the target (r > a) 

1 



Tl 



it2-h) = -1, 



and inside the target (r < a) 



—t2-(— + k]t 



Tl 



Tl 



Making a similar decoupling approximation as in two dimensions, leads finally to 

4. 



—lAt2--iti-t2) 

O T2 



These equations are solved inside and outside the target, using the following boundary conditions : 

-0, 



dr 



r—b 



dr 



dt 



2 

dr 



(A70) 
(A71) 
(A72) 

(A73) 

(A74) 

(A75) 

(A76) 
(A77) 
(A78) 



and the condition that t^l^) should be finite. 
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Figure 42 Static mode in 3 dimensions. \n{tm) as a function of ln(r2) for different values of n, a and h/a. Comparison b etween 
simulations (symbols) , analytical expression (line) and its asymptotics for fe ^ a obtained in (|Loverdo et al.L l2009bl ) (small 



dots), and simple expression for b > a and a small (|A79|) (dashed line), n ~ r"*" ~ 0.74^ ^ (|A80|) (blue, +), n = 0.25^ ^ 
(red, o), n = 2.5w^ (green, □). V = 1, k = 1. 



b. Results 

The following approximate expression of the mean search time is found in the low density limit 



fo'(T2 + ri) /(l + fcri) 



This expression of t„i can be minimized for : 



(A79) 



opt 



(A80) 



opt 



and the minimum mean search time finally reads 



'1.2 



(A81) 



topt-±i^ ( /«fc24i/4 , 5I/4 



(A82) 



Data obtained by numerical simulations (figure are in good agreement with the analytical expression (IA79[) 
except for small T2 or small b, where a refined analytical expression can be obtained (see ( Loverdo et al.l , [2009bl )). In 
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particular, the position of the minimum is very well approximated, and the error on the value of the mean search time 
at the minimum is close to 10%. 

With the static detection mode, intermittence is always favorable and leads to a single optimal intermittent strategy. 
As in one and two dimensions, the optimal duration of the relocation phase does not depend on k, i.e. on the description 
of the detection phase. In addition, this optimal strategy does not depend on the typical distance between targets b. 

One can notice than for the static mode in the three cases studied (1, 2, and 3 dimensions), we have the relation : 

opt / opt 



T-^ = Y /(2fc). This relation between the optimal durations of the two phases is independent from the dimension, 
b. Diffusive mode 

We now study the case where the detection phase is modeled by a diffusive mode. During the detection phase, the 
searcher diffuses and detects the target as soon as their respective distance is less than a (see figure I?!]) . 

a. Equations 

Outside the target (r > a), one has : 

^^•Vi2,0,0 + — (ti -t2.e,0) = -1, (A83) 
DAti + —(^[ ddsinO [ d(/)t2,e>-ti) = -1, (A84) 



and inside the target {r < a) : 

1 







"i^- Vi2,e,0 ^2,9,0 = -1, (A85) 

T2 

h = 0. (A86) 
With <2 = 47 Jq dOsinO J^^ ^0^2, 9,0, one obtains outside the target (r > a) : 

D A + ^ (t°"* - = -1. (A87) 

The decoupling approximation described in previous sections then yields outside the target : 

At™* + — (ir*-ir*) = -1, (A88) 

3 T2 

and inside the target (r < a) : 

A t^"* - = -1. (A89) 

3 T2 

These equations are completed by the following boundary conditions : 



dr 



= 0, (A90) 

r—h 



dt"^* 

dr 



tT\a)=iT\a), (A91) 

(A92) 



dr 
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b. Results in the general case 

Through standard but lengthy calculations the above system can be solved and leads to an analytical approximation 
of tm (see ( Loverdo et al.l , [2009bf ). In the regime h':^ a, and 6y^(r^I15)~M-^(T^w)~^ ^ 1; one obtains : 



KiK^TiDa |^tanh(K2a) + ^ tanh(K2a) 



with ^1 = ^;f;^^^ and «2 = It can be shown that i,„ only weakly depends on ti, which indicates that 

this variable will be less important than T2 in the minimization of the search time. The relevant order of magnitude 
for T°^* can be evaluated by comparing the typical diffusion length Ldiff = V 6Dt and the typical ballistic length 
Lbai = Vt. An estimate of the optimal time r°^* can be given by the time scale for which those lengths are of same 
order, which gives : 

« n 



In turn, the minimization of t„, leads to 



opt 



(A95) 



with X solution of : 

2 tanh(x) - 2a: + x tanh(a:)^ = 0. (A96) 



a 



This finally yields : 



Importantly this approximate expression is very close to the expression obtained for the static mode [t^ ~ \I Zv 

1.095p^) (jA81|) . and there is no dependence with the typical distance between targets h. The simplified expression of 
the minimal can then be obtained as: 



r^' ~ 1.078-. (A97) 

opt _ /6a 



and the gain reads: 



C* = {x - tanh(a:))-^ ^ 2.18^, (A98) 



gazn=^ ^0.15^, (A99) 

tni ^ 



where the search time without intermittence is given by tdiff ^ b^liiDa) for 6 > a (|Redneij . l200ll ) 



There is a range of parameters for which intermittence is favorable, as indicated by (see figure Both the 

analytical expression for tj^* in the regime without intermittence and with intermittence (jA98P scale as . However, 
the dependence on a is different. In the diffusive regime, tm oc a~^, whereas in the intermittent regime tm oc a~^. 
This enables us to define a critical Cn, such that when a > a^, intermittence is favorable: an — 6.5-^ is the value for 



which the gain (|A99|) is 1. 



c. Ballistic mode 



We now discuss the last case, where the detection phase 1 is modeled by a ballistic mode (see figure [2T|) . Since 
an explicit analytical determination of the search time seems out of reach, a numerical exploration of the parameter 
space is needed to identify the regimes where the search time can be minimized. Approximation scheme s are then 



developed in each regime to obtain analytical expressions (more details are given in (jLoverdo et al l l2009bl l 



a. Numerical study 

The numerical analysis performed in ( Loverdo et al.l[2009b[ ) puts forward two strategies minimizing the search time. 



depending on a critical value vf to be determined later on 

• when vi > vf, t°^* — > oo and Tj^* — > 0. In this regime intermittence is not favorable. 

• when vi < vf, t°^* — ^ 0, and r2^* is finite. In this regime the optimal strategy is intermittent. 
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Figure 43 Diffusive mode in 3 di mensions, tm /fdis as a function of r2 for different values of the ratio b/a (logaritlimic scale). 
The full analytical form given in (|Loverdo et al.l l2009bl ) (plain lines) is plotted against the simplified expression (|X93)) (dashed 
lines), the simplified expression with ti — (dotted line), and numerical simulations (symbols) for the following values of the 
parameters (arbitrary units): a — 1 (green, □), a — 5 (blue, *), a = 7 (purple, o), a — 10 (red, +), a = 14 (brown, x), a = 20 
(orange, o). n = 6 everywhere except for the small dots, V = 1, D = 1. tm/tdis presents a minimum only for a > ac — 4. 



b. Regime without intermittence (single state ballistic searcher) ; r2 — >■ 

Following the same argument as in two dimensions, without intermittence the best strategy is obtai ned in the limit 
Ti — » o o in order to minimize oversampling of the search space. Following the derivation of ()A64p (see ( Loverdo et al.l 
l2009bl ) for details), it is found that the search time reads : 

thai = (AlOO) 



c. Regime with intermittence 

In the regime of favorable intermittence, the numerical study suggests that the best strategy is realized for ti — > 0. 
In this regime ti — > 0, the phase 1 can be well approximated by a diffusion with effective diffusion coefficient 
DcS ~ vfTi/3 . The analytical expression tm derived in (jA93p can then be used, and yields for ti = Oand a : 

then one finds straightforwardly that Tj^* = where x is solution of xtanh(a;)^ + 2 tanh(x) — 2x = 0, that is 

X ~ 1.606. Using this optimal value of T2 in the expression of t™, one finally obtains : 

= ±^"^2.18-^. (A102) 




These expressions show a good agreement with numerical simulations (see ( Loverdo et al.l , [2009b( )) 



d. Discussion of the critical value Vi 
The gain is given by : 



5am -0.61-. (A103) 

tT VI 



As in two dimensions, it is trivial that vf < V, and the critical value can be defined as the value of vi such that 
gain = 1. This yields : 

~ 0.6V. (A104) 
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Importantly, Vi neither depends on b nor a. Simulations are in good agreement with this result, except for a small 
numerical shift. 



d. Conclusion in dimension three 

For the three possible modelings of the detection mode (static, diffusive and ballistic) in three dimensions, there 
is a regime where the optimal strategy is intermittent. Remarkably, and as was the case in one and two dimensions, 
the optimal time to spend in the fast non-reactive phase 2 is independent of the modeling of the detection mode and 
reads r2''* ~ l-l^- Additionally, while the mean first passage time to the target scales as Ir', the optimal values of 
the durations of the two phases do not depend on the target density a/b. 
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